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ABSTRACT 

In this paper, we demonstrate a new method for fitting galaxy profiles which makes use of 
the full multi-wavelength data provided by modem large optical-near-infrared imaging sur- 
veys. We present a new version of GALAPAGOS, which utilises a recently-developed multi- 
wavelength version of GALFIT, and enables the automated measurement of wavelength- 
dependent Sersic profile parameters for very large samples of galaxies. Our new technique is 
extensively tested to assess the reliability of both pieces of software, GALFIT and GALAPAGOS 
on both real ugrizYJHK imaging data from the GAMA survey and simulated data made to 
the same specifications. We find that fitting galaxy light profiles with multi-wavelength data 
increases the stability and accuracy of the measured parameters, and hence produces more 
complete and meaningful multi-wavelength photometry than has been available previously. 
The improvement is particularly significant for magnitudes in low S/N bands and for struc- 
tural parameters like half-light radius Tc and Sersic index n for which a prior is used by 
constraining these parameters to a polynomial as a function of wavelength. This allows the 
fitting routines to push the magnitude of galaxies for which sensible values can be derived 
to fainter limits. The technique utilises a smooth transition of galaxy parameters with wave- 
length, creating more physically meaningful transitions than single-band fitting and allows 
accurate interpolation between passbands, perfect for derivation of rest-frame values. 

Key words: methods: data analysis — techniques: image processing — galaxies: structure 
— galaxies: fundamental parameters 



1 INTRODUCTION 

Studies of galaxy formation and evolution rely on accurate esti- 
mates of physical galaxy properties, such as luminosity, mass, star 
formation history (SFH), size and morphology. Many of these prop- 
erties are obtained from imaging data, via measurements of mag- 
nitude and profile shape. Such measurements are nowadays rela- 



tively straightforward for individual objects. However, many anal- 
yses benefit from being applied to as large a sample as possible. 
Given the sizes of modem surveys, this typically means thousands, 
or even hundreds of thousands, of galaxies. For such large numbers 
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of measurements to be feasible, they must ideally be performed in a 
fully automated fashion, which significantly complicates mattersQ 

Most galaxy parameters may be obtained using a variety of 
methods. One common approach to measuring magnitudes is aper- 
ture photometry, in which one defines the extent of a galaxy in 
some manner and then sums all the flux within that area. Ideally 
one would choose an aperture large enough to contain effectively 
all the galaxy flux. However, one cannot simply use an arbitrar- 
ily large aperture, as that would introduce excessive noise from the 
sky and be more likely to be contaminated by flux from neighbour- 
ing objects. Typical methods employed to define photometric aper- 
tures therefore seek a reasonable compromise. As a result, aperture 
magnitudes necessarily miss a proportion of flux from the outer re- 
gions of each galaxy. Furthermore, the aperture defined for a given 
galaxy, and hence the amount of missing flux will vary between 
images, depending on spatial resolution, signal-to-noise (S/N) and 
the exact shape of the galaxy light distribution. Colour gradients 
within a galaxy also lead to the inferred extent of a galaxy vary- 
ing w ith observed wavelength (as recently reported bv lKelvin et alj 
120121 hereafter K12, and others). Further refinements include ap- 
plying convoluti ons to match the point spread functions (PSFs) of 
the images (e.g.. 'H ill et alj|2oTlh . and applying a minimal correc- 
tion by estimating t he flux that would be missed if the galaxy were 
a point source (e.g.. :White et alJl2005l ; lGraham et al.l2005l) . 

Varying definitions of photometric apertures can have a sig- 
nificant impact on the resulting science. Using a fixed surface- 
brightness threshold clearly misses more flux for objects that show 
less compact profiles. More sophisticated methods still suffer from 
biases, for example Petrosian magnitudes recover essentially all 
the flux for exp onential profiles, bu t miss ~ 20 to ~ 70 per 
cent of the flux i Graham et al.ll2005l) for de Vaucouleurs profiles 
jde Vaucouleursll948 f). The wavelength selected to define the aper- 
ture is also important. For example, disk galaxies are typically red- 
der in their centre, due to the presence of a bulge or dust. Defining 
the aperture in a red photometric band will therefore result in to- 
tal fluxes that are underestimated in bluer bands, and hence total 
colours that are systematically biased to redder values. 

For measuring galaxy sizes, one could employ methods simi- 
lar to those used to define photometric apertures. These approaches 
obviously suffer from many of the same issues described above, 
and generally do not provide a consistent, physically interpretable 
measure of galaxy size. A more meaningful alternative is to deter- 
mine the radius (or two axes of an ellipse) that contains a specified 
fraction of the total galaxy light. Common examples include the 
half-light radius, rso, and rgo, the radius which contains 90 per 
cent of the total galaxy light. Of course, these measurements de- 
pend critically on a reliable measurement of the total magnitude. 
Sizes derived using aperture magnitudes will suffer from system- 
atic biases with respect to galaxy profile shape, luminosity and dis- 
tance. There is again the issue of wavelength; sizes measured in 
blue bands will tend to reflect the extent of the young stellar pop- 
ulations, whereas red bands will more closely reflect the distribu- 
tion of stellar mass. Just as k-corrections are required to convert 
observed magnitudes to restframe values, similar corrections may 
be required to homogenise sizes when considering galaxies span- 
ning a range of redshifts. Finally, it is important to note that none 
of these size measurements are corrected for the effect of the PSF. 



^ Alternatively, for suitable tasks where automated tools are insufficient, 
one may inste ad resort to us ing large numbers of people, via citizen science 
methods, e.g. jLintott et aljf2008) . with their own set of complications. 



They will therefore be overestimated, particularly for intrinsically 
small or distant galaxies. 

A variety of automated proxies for morphology have been pro- 
posed, the simplest of which focus on the shape of the azimuthally- 
averaged surface-brightness profile. One widely used parameter 
is the concentration index, which is defined as the ratio of the 
radii conta ining two fractions of the total flux, e.g., C = rgo /rso 
(Strateva et al.l200l[) or C = 5 • log{rso/r2o) (as defined by CAS, 
iConselice<2003h . All of the biases which affect these size estimates 
will therefore result in biases in the concentration index and sim- 
ilar non-parametric profile measurements. In any case, even for 
well-resolved, bright galaxies, such simple proxies only give a very 
rough indication of true internal structure or morphology. 

MegaMorph is a project aimed at improving our ability to 
measure and understand the structure of galaxies. In particular, 
we endeavour to make optimal use of modem multi-wavelength 
imaging surveys. Using data from multiple bands simultaneously 
in the fitting process increases the signal-to-noise, without greatly 
increasing the number of free parameters. Importantly, combining 
multi-wavelength imaging provides information that is not avail- 
able to techniques which operate on only a single band. For ex- 
ample, this enables the fitting process to utilise the different wave- 
length dependence of each component to help separate their pro- 
files, and produces a more physically consistent models. We expect 
this to be particularly crucial when performing bulge-disk decom- 
positions. However, in this work, we first consider only fits using 
single-Sersic profiles. MegaMoiph and the software developed and 
utilized is further discussed in S ections [ 1 . 2 l and l2l 

1.1 Parametric methods 

To avoid many of the problems that empirical (aperture-based) 
methods suffer from, an increasingly popular approach to mea- 
suring galaxy properties involves fitting their surface brightness 
profiles with parametric models. This has a number of advan- 
tages: all the measurements are obtained in a consistent man- 
ner, varying PSFs can be easily accommodated, and the issue of 
missing flux is, at least partly, avoided. The price is the assump- 
tion of a parametric form for the two-dimensional surface bright- 
ness distribution; typically exponential profiles for galaxy disks, 
de Vaucouleurs profiles (de Vaucouleurs 19 48) for bulge s and el- 
lipticals, or more generally, Sersic profiles (ISersiclll968l) . A num- 
ber of software p ackages have been pro duced to perform such 
fits, e.g., GALHT (IPeng etalJl2002l,r2010l), GIM2d (ISimard|[l99^ 



'Simardetal.' '2002), BUDDA ('de Souza et al.' '2004'), 2DPH0T 
(LaBarberaetal. 2008) and GasPhot (Pignatelli et al. 2006). 

These tools can achieve good results, both when used man- 
ually to fit individual galaxies, or when applied to large surveys 
in a fully- or semi-automated fashion (at least w hen fitting sin- 
gle Sersic models to relatively bright galaxies, e.g.. iHaussler et all 
12007). However, while a user of this technique can successfully 
employ a complex combination of profiles whe n fitting individua l 
galaxies by hand (e.g., nearby NGC galaxies, Ivika et al]|2012h . 
automated model fitting in large surveys is considerably more chal- 
lenging. Many thousands of galaxies, each with their own indi- 
vidual complications, such as neighbouring objects and potentially 
varying sky level, PSF, imaging availability and profile complex- 
ity, must be dealt with in an automated fashion. Developing a fully 
automated code with sufficient complexity, accuracy, flexibility 
and speed to perform profile fitting in modem surveys is difficult. 
Neve r theless, a number of studies ( Schade et al. 1997 ; Lilly et al] 
1 19981; lAllen et alj|200^ ; ISimard et al.ll201 ll; iTasca & Whitell201 ll ; 
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I Kelvin et ail I2OI2I : iLackner & Gunr] l2012h have produced cata- 
logues of galaxy profile parameters for large samples, with recent 
local studies often bas ed on imaging pr ovided by the Sloan Digi- 
tal Sky Survey (SPSS. ' York et alfcOOOl) . Some of these works fo- 
cus on one-component Sersic models, although most also attempt 
bulge-disk decomposition. 

Many existing approaches to galaxy profile fitting are primar- 
ily designed to work with only a single image, i.e. only one pho- 
tometric band. If applied to a multi-band dataset, each band must 
be fit separately. One may choose to treat all bands equally, and 
allow the technique to fit a comp l etely independent model in each 
band (e.g.. lLa Barbera et al .l2010l:lKelvin et alj2012l) . This enables 
one to study wavelength-dependent structural variations, e.g., due 
to colour gradients. However, only a fraction of the data is used to 
constrain the profile in each fit. Furthermore, the resulting colours 
may not be physically meaningful, particularly in the case of mul- 
tiple components, due to unphysical variations in the structural pa- 
rameters (as demonstrated later). One might naively expect many 
model parameters to be identical in all bands, e.g., component cen- 
tres, axial ratios and position angles, or vary smoothly with wave- 
length, e.g., Sersic index and size. Alternatively, therefore, one may 
select one dominant band in which to fit an initial profile, and then 
fit this pro file to the other bands w hile holding various parameters 
fixed (e.g.. iLackner & Gunnll2012l) . With the profile fixed across 
all bands the resulting component colours should be more mean- 
ingful, but again only one band has been used to determine that 
profile, wasting data. Also consider that the smooth variation of pa- 
rameters with wavelength cannot be guaranteed which this method 
would assume to be the case. 

One approach to using all the available data to constrain the 
profile would be to simply sum all the images together and fit a 
model to the resulting image. Obviously, however, this does not 
allow colour information to be extracted. This profile could then 
be fit to the bands individually, with the structural parameters held 
fixed. Another solution is to fit a model to multiple images of the 
same object simultaneously. This is less common than fitting sin- 
gle images, but not a new idea. GIM2d ( Simardlll998 ) includes an 
option to fit two images with two bulge-l-disk models constrained 
to have the same structural parameters. Only the flux of each com- 
ponent is allowed to vary independently between the two models. 
This approach has been used to measure bulge and disk colours 
for over a million SDSS galaxies jSimard et~aLll201 ih . GIM2d also 
provides an ability to fit a stack of images with identical profiles, 
with only the centre of the model allowed to vary between images. 
However, we wish to (a) make use of an arbitrary number of multi- 
wavelength images, (b) constrain parameters to vary smoothly as 
a function of wavelength, being neither completely fixed or free, 
and (c) fit a variety of models, not just bulge-l-disk. We would also 
prefer to fit neighb ouring galaxies wher e appropriate, rather than 
relying on masking jHaussler et al.ll2007l hereafter H07). 

To understand the desire for model parameters which vary 
smoothly with wavelength, consider the example of fitting a single- 
Sersic model to a nomal disk galaxy, comprising a blue, ex- 
ponential disk and a red, de Vaucouleurs bulge. In bluer bands, 
the disk will be dominant, and hence the profile is best repre- 
sented by a low Sersic index, while in redder bands the bulge 
would become more dominant, resulting in a smoothly increas- 
ing Sersic index with observed wavelength (K12, Fig. 21). Addi- 
tionally, gradien ts in the stellar popul ations within spheroids (e.g.. 



La Barbera &d e Carvalho 2009; S uh et al il2010l) and discs (e.g. . 
Bell & de JongI I200Q. , Mac Arthur et~l.ll2004 iTortora et al.ll2010l : 



Gonzalez-Perez et al.ll201lh and centrallv concentrated dust atten- 



uation toriver et al.' '2007*; 'Maste rs et al.ll2010h can produce similar 
effects CMollenhoff et aL2006;.Pastrav et alj20T^ . Fixing the pro- 
file shape as a function of wavelength would therefore give bad fits 
in some bands. Allowing it to vary freely will often result in the 
parameters varying wildly with wavelength as the fit uses its in- 
creased freedom to fit the image noise. This will also significantly 
increase the number of parameters to be fit. 

We therefore propose that the preferred solution is to fit a 
full wavelength-dependent model to an arbitrary set of multi-band 
data, simultaneously. This approach can use all the available data 
to define the profile, while enabling the measurement of physically 
meaningful component colours and colour gradients. The form of 
the model parameters as a function of wavelength can be chosen to 
optimally balance consistency and flexibility. This approach should 
also improve the number of galaxies for which a full set of robust 
photometry can be determined, as the bands can 'help each other 
out'. For example, in a simultaneous multi-band fit, low S/N bands 
would not contribute much to defining structural parameters, but 
would benefit from the constraints on these from higher S/N bands, 
resulting in robust measurements of the flux in each band. 

1.2 The purpose of this paper 

In MegaMorph, we have developed a combination of tools in order 
to test our expectations regarding the benefits of multi-band para- 
metric measurements. This software is briefly described in Section 
[2] Details of the implementation, together with examples illustrat- 
ing the advantage of this approach, appear in Bamford et al. (2012, 
in prep; hereafter Paper I). Our technique is designed to be highly 
flexible, but for consistency we adopt a standard configuration for 
most of the work in this paper. Our choices are explained in Sec- 
tion[3] 

This paper is accompanied by another paper (Vika et al., in 
prep; Paper II), which applies our technique to a sample of 168 
nearby galaxies that have been artificially redshifted in order to 
assess its performance in fitting individual, realistic galaxies. The 
present paper complements that study by demonstrating the appli- 
cation of our technique to large surveys in an automated fashion, 
and with greater statistical power. Using both real and artificial im- 
ages, we will demonstrate how (and why) using multi-band fitting 
has advantages over single-band fitting, in tenns of stability, im- 
proved accuracy and increased sample sizes, especially for the low 
S/N bands of a survey. 

Most galaxies comprise multiple structural components, pri- 
marily a bulge and a disk. The most physically meaningful pa- 
rameters should therefore be obtained by fitting multi-component 
models. However, fitting such models is challenging, particularly 
on noisy, low-resolution, single-band imaging, as the parameters 
of the multiple components can be highly degenerate. Using multi- 
band data to constrain the fit significantly alleviates this problem 
(see Paper I), as the different wavelength dependencies of the in- 
dividual components (i.e. their colours) provides valuable informa- 
tion, which is not present in single-band fitting. 

Ultimately we aim to decompose galaxies into physically 
meaningful structures, and measure reliable properties for each 
component. However, simpler single-component fits still provide 
a great deal of useful information, and are less challenging to per- 
form. We will explore multi-component fits in future papers, but 
as a first step in demonstrating the advantages of multi-wavelength 
profile fitting, in our present work we focus on fitting single Sersic 
profiles to each object. 

This paper is structured in the following way: In § |2] we in- 
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troduce the idea of multi-band fitting, including a brief technical 
description on how this is carried out and what changes have been 
applied to both GALFIT (see § 12. Il l and GALAPAGOS (see § |2.2| l. 
§[3]explains the setup of both codes used throughout this paper. In 
§|4]we show tests from applying this software to real GAMA data 
(Galaxy And Mass Assembly) and compare the values to show in 
how much multi-band fitting improves the fitting results both on 
individual galaxies and on the galaxy population as a whole. § |5] 
carries out similar tests, but uses simulated images, e.g. galaxies 
whose true intrinsic values are known. This comparison, while not 
containing any real physical meaning about galaxy populations, al- 
lows to show the improvement by using multi-band fitting in more 
detail. §|6]takes other considerations than fitting accuracy into ac- 
count, e.g. fitting time and disk-space required. Finally, § |7] both 
as a sanity check and to further show improvements enabled by 
the new technique presented in this paper, discusses the colour- 
magnitude diagram of galaxies. This chapter is aimed to be a mo- 
tivation for users to apply the software developed, tested and pre- 
sented in MegaMorph papers in order to improve their scientific 
results. 



2 MULTI-WAVELENGTH PROFILE FITTING 

In order to evaluate the advantages of fitting wavelength-dependent 
models to multi-band data, we have implemented software to per- 
form such fits. For the sake of efficiency and reliability, we chose 
not to re-implement all the functions required for a profile fitting 
code ourselves. Instead we elected to build upon existing, well- 
tested software and make only those changes necessaiy to enable 
multi-wavelength fitting. However, in the course of modifying the 
software, we have also added additional features where required or 
deemed convenient, and generally improved the efficiency of the 

code where possible. 

We selected Galapagos jBarden et ar]|2012l) and GALFIT3 
jPeng et alJlioToh as the starting point for our development, due 
to their reputation for reliability, flexibility and speed, as well as 
the extensive experience of members of our team in using these 
software tools (H07). GALFIT performs the fit for each target image 
while taking the image PSF into account; GALAPAGOS, after initial 
preparation of the data, takes care of everything else required to run 
GALFIT in an automated manner on a large surve y, including book- 
keeping, object detection (using SEXTRACTOR: lBertin & AmoutsI 
[l99 6), cutting images of each target, masking, determination of the 
sky level, estimation of initial parameters, writing setup files and 
load-balancing. In this section we briefly describe our choices in 
implementing multi-wavelength fitting, and the modifications we 
have made to the standard versions of these codes. 



2.1 GALFITM 

We have adapted GALFIT3 jPeng et al]l20I(]|) for the requirements 
of this project, with permission of the original developer, C. Peng. 
To differentiate our modified version from the standard release we 
refer to it as GALFITM. For reference, all the work in this paper uses 
GALFITM version 0. 1 .2. 1 . The code will be publicly released in the 
near future. GALFITM-0. 1.2. 1 is based on GALFIT version 3.0.2, al- 
though the additions in GALFIT-3.0.4 (the latest standard version) 
will be incorporated in GALFITM before public release. Develop- 
ment is continuing, primarily to improve ease-of-use and incorpo- 
rate the additional features mentioned above. However, the general 



performance of the technique is expected to remain as presented in 
this paper 

GALFIT constructs model images by summing one or more 
components, which potentially include a sky background (with op- 
tional gradient), elliptical Sersic functions, point sources and a va- 
riety of other profiles. GALFIT fits the parameters of its model to 
the input data (weighted by an ertor map, which may be provided 
or internally created) by employing the widely-used Levenberg- 
Marquardt (LM) algorithm to minimise the weighted sum of the 
square residuals (x^)- In addition to the model image itself, GAL- 
FIT calculates the derivative of the model image with respect to 
each free parameter, as required for the LM algorithm. The model 
and all of its derivatives are convolved with the provided PSF for 
comparison with the in put image. The reader is advised to consult 
IPeng et al.l ( l2002l.l2010i) for a detailed description of GALFIT. 

The standard version of GALFIT3 accepts only a single input 
image with which to constrain the model fit. It was therefore nec- 
essary to make fairly substantial modifications to enable the use of 
multi-band data. However, most of the original code and its struc- 
ture is maintained, and we intend our modified version to be back- 
ward compatible when used with single-band data. In this subsec- 
tion, we briefly describe the significant changes. For full details we 
refer the reader to Paper I. 

2.1.1 Wavelength-dependent model parameters 

In order for GALFITM to be able to fit multi-band data, we replaced 
every galaxy model parameter with a wavelength-dependent func- 
tion, 

M{x,y; pi, . . . ,p„) M[x,y-pi{X; 171,1, . . . ,gi,„J,. . . , (I) 

Pn(A; q„,i, . . . ,(7n,m„)] , 

where M{x,y; ■) is the model (describing the surface brightness 
as a function of pixel coordinate, before PSF convolution), the pi 
are the n original parameters of the GALFIT3 model, and each pi is 
some function, with m,i parameters g^j, which describes the varia- 
tion of the model parameter i versus wavelength, A. Whereas stan- 
dard GALFIT fits the Pi, in GALFITM the parameters of the fit are 
the set of qi^j. 

In the case of a standard Sersic profile used in this paper, these 
parameters are position, magnitude, half-light radius, Sersic in- 
dex, axis ratio and position angle. However, the approach is imple- 
mented in a general fashion and works for any of the model func- 
tions provided by GALFIT3. The choice of function is somewhat 
arbitrary, although various properties are clearly desirable, includ- 
ing a straightforward way of selecting the function's flexibility and 
independence of the function parameters. We chose to use a series 
of Chebyshev polynomials (of the first kind; lAbramowitz & Steguj 
IT965.) , Tj (z), for all of the functions: 

m j 

?5.(A; te.i}) = ^*.,r,[z(A)]. (2) 

j=o 

The Chebyshev polynomials are restricted to the domain 2 G 
[—1, +1], and hence the wavelength range of the input bands is lin- 
early mapped on to that interval. The Chebyshev polynomial Tj (z) 
is of order j, i.e. T2{z) is a quadratic function of z. Cheybshev 
polynomials of the first kind are defined by the recurrence relation 

To(z) = 1 
Ti(z) = z 

T„+x(z) = 2zT4z) - r„-i(z) . (3) 
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The fit parameters, qi,j , are therefore the Chebyshev coefficients of 
the series for parameter i. 

The flexibility of this function may be varied by selecting the 
maximum order of each series, rrii, i.e. limiting to zeroth-order im- 
plies that a parameter must be constant with wavelength, second- 
order allows quadratic dependence with wavelength, while choos- 
ing the order as one less than the number of bands gives the function 
freedom to interpolate the data precisely. 

Chebyshev polynomials fom an orthogonal basis set, but full 
orthogonality only occurs when the function is constrained at the 
corresponding set of Chebyshev nodes. This is possible when ap- 
proximating a smooth function, but in our case we are not free to 
choose the position of the constraints, they are set by the locations 
of the available set of photometric bands on the selected wave- 
length scale. Nevertheless, the independence of individual Cheby- 
shev polynomials, even if only partial, is expected to limit the de- 
generacies between parameters, and hence aid the stability of the 
fitting process. 

We stress that the purpose of the functions, pi (A), is to connect 
the parameter values in the different bands with a user-specified de- 
gree of smoothness. For example, we might choose the Sersic index 
to vaiy quadratically, position angle constant, and magnitude to be 
completely free as a function of wavelengtf0. The functions them- 
selves are not intended to be physically meaningful, although they 
may be used to approximate parameter values at wavelengths be- 
tween the observed bands, e.g. to determine restframe parameters. 
In this paper we use the Chebyshev polynomials as a function of 
wavelength. However, the variable used in the polynomials need 
not be true wavelength. Frequency, the logarithm of wavelength, or 
a variety of other variables may be suitable (see Paper I for a more 
thorough discussion). For the purpose of this paper, we have chosen 
to use linear scaling with wavelength. 

All of the free parameters of the model, the set of qtij, are fit 
to all the multi-band data simultaneously by minimising a single 
quantity, defined as: 

^2 ^ [du,v,u, - M{xu,yv\\u,,qi.j)f 

u,v,w "u,v,w 

where it and v index the pixels, at positions Xu and yu , in image , 
with wavelength A^,, and uncertainty image a^,. We write Eqn.|4] 
in this way to emphasize that the data comprises a set of discrete 
measurements, while the model M is, in principle, a continuous 
function, evaluated at the position and wavelength of the data in 
order to compute x^. Further technical details will be presented in 
Paper I. 

2.1.2 Parameter constraints 

One side effect of our multi-wavelength modifications is that the 
approach taken to constrain model parameters in GALFIT3 required 
revision. These constraints take two forms: hardcoded limits (such 
as ensuring that sizes cannot become negative) and user specifiable 
limits, but both are treated similarly. Constraints are useful to guide 
the fitting process, by eliminating regions of parameter space which 
are ruled out by other considerations. They can therefore improve 
the efficiency of the early stages of the fitting process. However, 
if the fitting process repeatedly encounters constraints, this is an 
indication that a good model fit to the data cannot be achieved. 

To be exact, we use 3543A, 4770A, 623lA, 7625A, 9134A, 10305A. 
12483A, 16313A, 22010A for ugrizYJHK-hand, respectively. 



Appropriate handling of constraints is particularly important 
when using GALFITM to fit multiple objects simultaneously (as is 
conmion with GALAPAGOS). When considering a single galaxy 
(possibly with multiple components), if the fitting process ends 
with a parameter very close to a constraint boundary, it is reason- 
able to discard the resulting fit from subsequent analysis. How- 
ever, in the case of a target galaxy with one or more neighbours, 
we would not want difficulties encountered in obtaining an uncon- 
strained fit for a neighbouring object to negatively impact the fit to 
the primary target, or result in a potentially good fit to the primary 
target being discarded. 

In GALFIT3, the physical parameters, pi, feature directly in the 
fitting algorithm. Constraining these physical parameters to lie on 
specified intervals can therefore be achieved in a straightforward 
manner. At each iteration of the fitting process, the LM algorithm 
proposes a step for each parameter. If that step would violate a con- 
straint, GALFIT3 typically resolves the conflict by simply setting 
the offending parameter to the value at the constraint boundary. 
For multi-band fits, however, constraints on the physical parame- 
ters may be violated in some bands but not others. There is also a 
complicated relationship between the physical parameter at a given 
wavelength, Pi(\), and the fit parameters, qi,j. An alternative ap- 
proach is therefore required. We briefly outline this here. For fur- 
ther details and discussion see Paper I. 

The LM algorithm interpolates between the Gauss-Newton 
(GN) algorithm and the method of gradient descent (GD), with the 
degree of interpolation controlled by a damping parameter, A. The 
GN algorithm will generally attempt to make relatively large steps, 
whereas GD is more conservative. Increasing A leads to dominance 
of GD over GN, and increasingly smaller steps. The LM algorithm 
includes a prescription for varying A to appropriately balance GN 
and GD as the fit progresses: if a proposed set of parameter steps 
successfully improves , then A is decreased by a factor, and the 
steps are accepted; otherwise it is increased by the same factor, and 
the steps are rejected. (This factor is 10 in the case of GALFIT.) 

In GALFITM, if a proposed set of steps in the fitting parame- 
ters Qi^j would violate a constraint on the standard parameter pi in 
any of the wavelength bands, then the steps are not performed for 
parameters g^,.. All other (unoffending) parameters are stepped as 
usual and a trial value of generated. This approach avoids the 
difficulty of determining how to limit the parameters g^,. so as to 
avoid the resulting Pi{\; {qi, }) from violating any constraints at 
the wavelengths of the input bands. 

For a moment assume that this is the only change. In that case 
a failure to improve x^ would result in a decrease in A, and re- 
jection of the entire proposed set of parameter steps. The steps in 
the next proposed set would be smaller, and more likely, though 
far from guaranteed, to avoid violating constraints. However, often 
the proposed steps (without those which would cause violated con- 
straints) will result in an improved x^ , an acceptance of those steps 
and an increase in A. The next iteration is therefore likely to pro- 
pose a step in the offending standard parameter that is similar to, 
or larger than, the previous. Tests have shown that this often results 
in that parameter remaining fixed for the entire duration of the fit, 
even though a small movement toward the constraint may lower x^ 
and result in a more acceptable model. 

To mitigate this issue, and encourage movement of constraint- 
violating parameters towards (but not beyond) constraint limits, we 
impose a schedule for A. This is designed to substantially increase 
A occasionally in the case of violated constraints, resulting in the 
next set of proposed steps being much smaller than the previous 
proposal, and so more likely to avoid overstepping the constraint 
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boundary. This adopted schedule was developed through trial and 
error, but appears to do a reasonable job of meeting our require- 
ments. 

The result is that constraints can be specified in GALFITM, 
on both individual fitting parameters, qi,j, and more usefully on 
the standard parameters, pi, at any, or all, of the input band wave- 
lengths. For example, in the fits in this paper the Sersic index, n, is 
constrained to lie on the interval 0.2 ^ n ^ 8 at all input wave- 
lengths. Details of all the constraints applied in the present work 
are given in Section|3] 

2.1.3 Other modifications 

While making the changes described above, we have attempted to 
retain backward compatibility as far as possible. While the setup is, 
of course, slightly more complicated for multi-band data, we ensure 
for any original GALFIT3 start file (for single-band data) to work 
unaltered with GALFITM. The additional multi-wavelength features 
are simple to enable, and any user already familiar with GALFIT 
should have no difficulties using GALFITM. 

The output format of GALFITM is also slightly modified. In 
addition to the image, model and residual (which of course GAL- 
FITM provides for each band), for convenience it also stores the 
PSFs used and provides all the fitting information, including setup 
details and the full results, in FITS tables within the output file. In- 
formation is also still provided via header keywords for backward 
compatibility. 

Several minor fixes and efficiency improvements have also 
been made. For example, all variables are now stored as double 
precision. Besides being more accurate, this also provides a mod- 
est speed improvement on modeiTi 64-bit machines. Again, we refer 
interested readers to Paper I for full details of our GALFIT modifi- 
cations. 

As part of our MegaMorph project we are investigating several 
other modifications to GALFIT, including the incorporation of non- 
parametric components, and alternatives to the LM algorithm. The 
results of these investigations will be described in future papers. 

Throughout the development of GALFITM, we have compared 
its output on single-band data to that of GALFIT3, generally finding 
very close agreement. For the vast majority of galaxies, the results 
of GALFITM and GALFIT3 are identical. The greatest differences re- 
late to our modified implementation of constraints. In cases where 
constraints are encountered during the fit, this can result in some- 
what slower convergence, but this is often accompanied by a for- 
mally better fit, in terms of a slightly reduced compared with 
GALFIT3. 

2.2 GALAPAGOS 

We have adapted the current public version of the IDL script GALA- 
PAGOS (version 1.0) to support GALFITM, and hence utilise multi- 
wavelength data, in close collaboration with its original developer, 
and coauthor of this paper, M. Barden. In the process, we have 
implemented a number of improvements and additional features, 
some for efficiency or convenience, and others that were required 
by the nature of our chosen dataset. We refer to our new version 
as GALAPAGOS-2. Specifically, version 2.0.2 was used to produce 
most of the results shown in this paper. Section|7]uses version 2.0.3, 
but the two versions only differ in very minor details and should 
produce nearly identical fitting results. Changes became necessary 
in order to be able to target specific objects instead of every ob- 
ject that was detected, hence speeding up the analysis for Section 



7, where only a subset of the detected objects are considered. In 
all versions of the code, we have attempted to preserve backward 
compatibility in the case of single-band data. The code will be pub- 
licly released in the near future. In this subsection, we will briefly 
explain our modifications. 

Prior to this work, GALAPAGOS was designed for use with 
space-based imaging, specifically for surveys performed by the 
Hubble Space Telescope (HST). The HST PSF is very stable on 
the multi-drizzled images that are usually used for fitting purposes, 
both temporally and across the field of view. GALAPAGOS- 1 there- 
fore only used one PSF for the entire survey. For the MegaMorph 
project, we wished to apply GALAPAGOS to ground-based imaging. 
While it is possible to approximately homogenise the PSF across an 
entire survey, by applying appropriate convolutions, the resulting 
resolution must necessarily correspond to the worst-case and thus 
a great deal of spatial information would be discarded. Instead, we 
adapted GALAPAGOS to work with a spatially variable PSF. 

In principle, it would be desirable for GALAPAGOS to con- 
struct an estimated PSF for each target galaxy (as was done in K12). 
However, this would require providing GALAPAGOS with knowl- 
edge of the survey strategy. As each survey will typically adopt 
a different tiling strategy, PSF creation is not trivial to generalise, 
particularly given the importance of the PSF in correctly modelling 
galaxy profiles. Sophisticated software already exists to take a set of 
point sources and combine them to produce an accurate PSF (e.g., 
PSFE)13) and a user of our code should use those to pre-determine a 
set of suitable PSFs. For the GAM A survey, PSFs had already been 
determined by K12 and were used throughout this analysis. 

For generality, we implemented a selection of the PSF from a 
provided list of filenames and sky coordinates. GALAPAGOS sim- 
ply selects the closest PSF to each target position. As K12 generate 
PSFs at the position of each galaxy with rrir < 19.8, the PSFs 
will correspond exactly for these galaxies, and the sampling is suf- 
ficiently dense that, in the vast majority of cases, fainter galaxies 
will be well represented by their nearest PSF. The PSF selection is 
performed for each band individually. 

We ultimately aim for our technique to be applicable to the 
largest surveys available. During this proof-of-concept stage we 
are content to restrict ourselves to more modest datasets, but still 
wished to work with a single GAMA II region, of area ~ 60 deg'^ . 
We therefore improved the efficiency of the code in several places 
where it became apparent that, for large datasets such as ours, a 
major speedup was possible. Some of these changes may alter the 
outcome of the code very slightly (e.g., sub-scripts of GALAPAGOS 
now know only about the neighbouring frames at times, instead of 
the entire survey, potentially changing deblending decisions). How- 
ever, we think that in practice this will not produce any significant 
differences in the results, as we were always very conservative in 
our modifications. Overall, we were able to speed up the code by 
approximately a factor of four in terms of CPU time. A further sim- 
ple optimisation was made in the loop over all objects in the sur- 
vey. GALAPAGOS- 1 determines whether the next object in queue 
is close enough to be influenced by an object currently being pro- 
cessed, and, if so, waits for that object to finish before starting the 
next. Our new code simply continues with a different object further 
away, thus keeping more CPUs busy at any given time. 

Most importantly, of course, the code can now handle multi- 
wavelength datasets. In addition to defining all the bands, each 
comprising a set of images, which are to be used for fitting (using a 
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setup similar to GALAPAGOS- 1), tlie user must define an additional 
set of images on which SEXTRACTOR is run for object detection. 
Of course, the detection images could simply correspond to one of 
the fitting bands. In our case, we chose to use a co-added image 
of all bands for detection, in order to detect sources with extreme 
colours that would be missed if only one band were used for de- 
tection. In this way, we use the whole dataset in order to derive an 
object list. Even so, the fitting bands are not all quite equal: one 
must be defined as the primary band, on which all deblending and 
masking decisions are made. We chose the r-band in our dataset, 
because it typically has the highest S/N. Our construction of the 
multi-band detection image is very simple, and could admittedly 
be i mproved upon by w eighting the input images more carefully 
('e.g.. lSzalavetaD fT99i). Nevertheless, it suffices for the purposes 
of this paper. 

The output returned by GALAPAGOS-2 provides much more 
information than previously; namely most of the details from the 
FITS tables in the GALFITM output file. However, GALAPAGOS-2 
may continue to be used with GALFIT3, for single-band data, in 
which case it provides the same information as the previous ver- 
sion. 

Further planned changes include the implementation of bulge- 
disk decomposition, ideally together with model selection, so that 
GALAPAGOS itself is able to decide whether a single-profile fit 
or bulge-disk decomposition provides the better representation of 
the imaging for each object and hence the more useful set of pa- 
rameters. We also plan to adapt GALAPAGOS to run in supercom- 
puting environments, in order to achieve the speed necessary for 
larger samples and/or surveys. However, the general performance 
of single-Sersic profile fitting is expected to remain as presented in 
this paper. 



3 CHOICE OF MULTI- WAVELENGTH MODEL, INITL4L 
PARAMETER VALUES, AND CONSTRAINTS 

GALAPAGOS requires various choices to be made regarding its op- 
eration and the setup information it provides to GALFIT. Our gen- 
eralisation to multi-band data adds a number of additional options. 
This section describes the choices we have made for the analysis 
described in this paper. 

If GALFITM is used to fit a single-component profile with 
structural parameters (i.e., all except magnitude) that are constant 
with wavelength, this is mostly equivalent to using GALFIT3 to first 
fit a single co-added image to obtain these parameters, and then 
measuring each magnitude by fitting each band with this fixed pro- 
file. When using multiple-component profiles, this is no longer true, 
and retaining the multi-band information throughout the fit leads 
to more accurate and reliable measurements (see Paper I). In the 
case of single-component fits, the greatest advantage of our new 
fitting technique comes from allowing profile structural parameters 
to vaiy systematically with wavelength; the use of co-added im- 
ages would lose this information. We will quantify the benefit of 
our multi-band fitting approach in Section |5] An important set of 
choices are therefore the degree of wavelength dependence we al- 
low for each parameter. 

It is critical that we obtain an accurate magnitude for each 
band, and hence colours. The magnitudes in each band are clearly 
correlated, and the variety of possible galaxy spectral energy dis- 
tributions (SED) are very well known. These cannot be reproduced 
with a low-order polynomial, and so we must ensure that sufficient 
freedom is given to the magnitudes such that they are accurately 



recovered. Full freedom is implied by using a polynomial with as 
many coefficients as data points, in which case the function is capa- 
ble of perfectly interpolating the data. To describe the wavelength 
dependence of magnitude, we therefore use an 8th-order polyno- 
mial (with 9 free coefficients, equal to the number of bands in our 
dataset). Note that the use of high-order interpolating polynomials 
is afflicted by Runge's phenomenon, whereby the function oscil- 
lates excessively between data points, particularly at the edges of 
the considered interval. The entire function itself therefore does not 
well-represent the galaxy SED and so cannot be used to estimate 
magnitudes at wavelengths other than those for which there is con- 
straining data (see Section lS^ . However, the magnitudes obtained 
for each band in the dataset remain reliable. If required, e.g., for 
k-correcting magnitudes, the resulting magnitudes may be interpo- 
lated using different codes, based on r ealistic SED templates (e.g., 
KCORRECT; lBlanton & RoweislboOTi) . Finally, note that while we 
allow full freedom for magnitudes in this paper, in further work it 
may be appropriate to consider using slightly lower order polyno- 
mials to reduce this issue, while still retaining sufficient flexibility 
to recover accurate values (see Paper I). 

Equally important to accurate magnitudes is the determination 
of physically meaningful structural parameters of the galaxies, e.g. 
half-light radius and Sersic index. We chose to allow profile half- 
light radius and Sersic index to vary with wavelength quadratically 
(i.e., second-order polynomials, with three coefficients). This was 
decided after examining the wavelength dependence of these quan- 
tities from single-band fits to bright galaxies. For most of these a 
linear function was sufficient to model the trend, but in some cases 
a mild curvature was seen. We therefore elected to fit a polynomial 
one order higher than linear, in order to examine this effect. 

In our simulations (described in Section |5j, we create simu- 
lated galaxies for which half-light radius and Sersic index vary ac- 
cording to known second-order functions. When fitting these simu- 
lations, we allow these parameters to vary with third-order, in order 
to investigate our ability to recover the correct higher-order coeffi- 
cients. 

As the images for each wavelength band are accurately regis- 
tered (although read Kelvin et al. 2012 and §4.1), the centre of the 
profile is constant with wavelength (although the centre is allowed 
to vary during the fit; we do not constrain the position to that given 
by SEXTRACTOR). Similarly, for position angle and axis ratio we 
also choose to fit constant values, with no wavelength dependence. 
While this ignores variations that might be expected for typical spi- 
ral galaxies (red bulges are round, blue disks appear elongated), this 
seemed to be a reasonable approximation for our purposes in this 
paper. 

The sky values for each band are pre-determined by GALA- 
PAGOS and held fixed during the fifl. We have shown in H07 that 
this is the most reliable approach for single-band fits, and there is 
no obvious reason why this should not also be the case when using 
multi-band data. 

Largely following H07, we adopt the following constraints on 
the parameters in GALFITM. In the case of multi-band fitting, these 
constraints apply to the parameter values for all bands (but not the 
entire polynomial). 

* At this point, it should be mentioned, that the GAMA data shows im- 
perfect flat-fielding of the provided Swarped images due to the use of lai'ge 
filters. These are chosen to avoid removing too much real local stracture 
but lead to sky backgrounds not being accurately measured around small 
objects. As this effect would be present in both single and multi-band fit- 
ting, we ignore this effect in this paper. 
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Position (x, y): These are simply constrained to lie within the im- 
age cut-out for this object. In practice, this constraint is rarely en- 
countered during the fit, but is retained to prevent the centre running 
out of the image in the case of a nearby bright source. 

Magnitude (m): —5 ^ mflt — minput ^ 5, where is de- 

rived by adding an empirically estimated offset to the Mag_Best 
derived by S EXTRACTOR during object detection. Additionally, we 
use ^ m ^ 40, to ensure sensible values; in practice this con- 
straint is rarely hit. 

Size (half-light radius; ro): 0.3 ^ re ^5 400 pixels. This main- 
tains values in a physically meaningful range and prevents the code 
from fitting very small sizes, where, due to oversampling issues, 
the fitting iterations become very slow. Pixel sizes in the data used 
are 0.339 arcsec/pixel, hence we constrain the half-light radii to be 
rc > 0.1 arcsec. For reference, it should be noted that re > 0.1 
arcsec corresponds to re ~ 0.3 kpc at a typical objects redshift of 
z ~ 0.2 in the GAMA survey. 

Sersic index (n): 0.2 ^ n ^5 8. Fits with values outside these 
ranges rarely represent good models of a target galaxy. The upper 
value of 8 is a conservative choice as objects with higher Sersic in- 
dices are rarely seen and, from earlier visual inspection, are usually 
associated with spurious galaxy fits or cases where the target object 
is a star It should be sta ted that some lumin ous elliptical galaxies 
with n > 8 do exist (e.g. lGraham et al .l2005h . hence this constraint 
will be removed and loosened to higher values in future works. 

Axis ratio (g): 0.0001 ^ q ^ 1. Again, this ensures the fit value 
is physically meaningful, but is mostly superfluous, as GALFIT in- 
cludes a hardcoded constraint on ^ g ^ 1. The main reason 
for this constraint being applied is that GALFIT, when very small 
values are reached becomes very slow due to oversampling issues. 

Position angle (61): — ISOdeg < 6 < 180 deg. This constraint 
is hardcoded in GALFIT. Following the same definition as GALFIT, 
position angle is defined as a major axis positioned vertically is 
deg (nominally north if rotated to the standard orientation) and 
increases counterclockwise (nominally toward the east). 

These constraints are implemented in order to improve and 
speed up the fitting process when fitting galaxies. For stars, the sit- 
uation is slightly different. An unsaturated star should technically 
return a point source when PSF-correction is used during the fit, 
i.e., re = 0. However, for practical reasons, enforced by the con- 
straints specified above, the fit is not allowed to do this. Instead, 
the fit results usually end up on one of these fitting constraints (typ- 
ically re = 0.3[pix] and n — 8) and thus removes the star from 
the image in a slightly non-optimal way. However, when this model 
is subtracted from the image after PSF convolution, we generally 
find these constraints produce a reasonable residual image, thus not 
significantly influencing the galaxy fits. For saturated stars, fitting a 
Sersic fitting, while masking out the saturated part of the profile, is 
more suitable to remove the wings of the profiles, more important 
in order to improve the fit of the neighbouring galaxies. 

We use these constraints to remove both stars and galaxies 
with bad fits, from our catalogue, by identifying objects with fit re- 
sults lying on one or more constraint boundaries. As just explained, 
the vast majority of stars result in values on one of the above con- 
straints. This also occurs for galaxies when the object in question 
cannot sensibly be fit with a Sersic profile, and hence any returned 
values should not be used in a scientific context. Please read Sec- 
tion |4j2] for more details about the cleaning of the fitting results 
catalogue. 

One obvious potential improvement would be to fit stars with 
PSF profiles (as has been done in K12) instead of Sersic profiles. 



GALFIT generally does allow this, but there are two main reasons 
why we have chosen not to do this. Firstly, such a procedure would 
require a reliable galaxy/star classifier, to make the decision of 
which profile to use for each object. While this feature is desir- 
able, it is not straightforward to implement, and was not deemed to 
be high-priority, given the low impact we expect it to have on our 
results. Secondly, there are bright stars present in the images, which 
possess significant flux at radii beyond the size of the PSF im- 
ages we use. This is especially true for very bright, highly- saturated 
stars, for which the profile core strongly deviates from the PSF due 
to saturation effects. Furthermore, the wings of bright stars often 
vary, and are therefore not well represented by an averaged PSF. In 
these worst cases, it was found that fitting a Sersic profile, rather 
than a PSF, results in much better residual images, particularly at 
large radii, as it has more flexibility to mimic and remove the outer 
wings of the PSF profile. In our data setup, we additionally use a 
masking scheme to identify saturated areas in the image. These ar- 
eas are consequently (a) smoothed in the S EXTRACTOR detection 
image, such that star images with the internal structure typical of 
saturated sources in SDSS are detected as one object only; and (b) 
masked in the images used by GALFIT, hence sufficiently removing 
the wings of the stars and masking out the innermost areas, so a fit 
to a neighbouring source should not be significantly influenced by 
either of these areas. 



4 APPLICATION TO REAL IMAGING 

We evaluate our multi-band galaxy profile fitting technique on both 
real and simulated datasets. Our real data obviously have the ad- 
vantage of showing actual galaxies, and we can compare to results 
from other studies. However, there is no definitive 'truth' to which 
we can compare. Simulated data, on the other hand, are idealised 
and do not capture all the subtleties of real data, and of course can- 
not be used to study the real universe. However, they do provide 
us with a way of testing the results of our method against known 
values. We first present our work using the real dataset. 

In this section, we describe the data, and show the results of 
extensive tests and comparisons between single-band and multi- 
band fitting techniques, using an otherwise identical code. 

4.1 Data 

For the purpose of these tests, we have chosen to u se imaging data 
provided by the GAMA survey l iDriveretal.l201 ij), as it comprises 
one of the largest multi-wavelength datasets currently available, in 
terms of both area and wavelength. GAMA is focussed around a 
redshift survey, but, crucially, this is supplemented by a highly ho- 
mogeneous and complete set of multi-wavelength data, spanning 
from the far-UV to radio, making it a superb tool for studying the 
local universe. We plan to use the galaxy profile fits from this paper 
and future work to perform a variety of scientific studies for which 
the wealth of the GAMA dataset is extremely valuable. 

The present GAMA survey (phase I) covers an area of ~ 144 
deg'^, of which half is in three ~ 48 deg^ equatorial fields, with high 
spectroscopic completeness to a depth of r < 19.4. Current efforts 
(GAMA II) are focussed on two additional, more southerly, fields 
as well as expanding the equatorial fields to ~ 60 deg^ , increasing 
the depth of all fields to r < 19.8 and the total survey area to 
~ 290 deg^. As we only require a relatively modest sample size 
for the purposes of this paper (~ 10* galaxies), and because of the 
computing time required to fit large samples, in this paper we limit 
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ourselves to one region of GAM A, the equatorial field at 9h R.A., 
and often to only a subregion of that. In future work we will expand 
our analysis to the full GAMA II survey. 

GAMA has prepa red its imaging data in a very convenient 
foiTn for our purposes dHill et al.ll201 ll) . These data include five- 
band optical (ugriz) imaging from SDSS plus four-band near- 
infrared (YJHK) imaging from the Large Area Survey (LAS) 
component of the U KIRT Infrared Deep Sky Survey (UKIDSS; 
[Lawrence et al.l 12007 ). All of these bands have a depth and reso- 
lution amenable to Sersic-profile fitting. Importantly, the images 
for all nine bands ha ve been 'micro-re gistered' onto the same pixel 
grid, using SWaip teertin et al. This procedure also ho- 

mogenises the photometric zeropoints and sky background in order 
to p roduce an artefact-fr ee mosaic. 

iKelvin et alj j20I2h has presented the results of single-band 
Sersic-profile fitting for all GAMA spectroscopic targets (r < 
19.8) in the GAMA I fields. This catalogue has provided a very 
useful comparison dataset, enabling us to carry out initial tests and 
aiding the early development of our code. We have used our own 
software to perform single-band fits, to ensure a fair comparison 
with our multi-band results, and so fit parameters from K12 are not 
shown in this paper other than in Fig. [T] A more complete compar- 
ison between the two different single-band codes used here and in 
K12 is beyond the scope of this paper and requires a detailed dis- 
cussion about minor details in the two codes. However, we have 
made detailed comparisons between our single-band results and 
those of K12, finding generally excellent agreement. The PSFs used 
in this paper are those obtained by K12 with GALAPAGOS' choos- 
ing the closest PSF to the targeted object for the fit. 

The imaging data used in this analysis will be made public in 
GAMA data release 2. Compared to the already public DRl im- 
ages, these images cover a larger area with a slightly different pixel 
scale to better match forthcoming GAMA datasets. The only mod- 
ifications we make to these imaging data is the addition of a back- 
ground pedestal so that GALFIT can construct correct sigma images. 
For each band, we simply use a typical background value from the 
original SDSS imagingfj. We also cut the images into overlapping 
tiles, which makes the data easier to handle in GALAPAGOS. 

A second set of images is used within GAMA to carry out the 
aperture photometry used in Section |7] In order to derive aperture 
matched photometry, all images have been blurred to a common 
PSF size of 2 arcseconds FWHM, creating a homogeneous dataset 
and avoiding artificial colour biases that would be introduced by 
the individual seeings in the 9 observed bands. 

4.2 Profile fitting 

We used our multi-band version of GALAPAGOS (itself utilising 
GALFITM), to create several catalogues for all objects detected in a 
small subregion (~ 1.3° x 1.3°) of the 9-hour GAMA field (G09) 
using the same versions of our codes to ensure maximum com- 
patibility of the results. In all runs on the data, we have used the 
same setup as closely as possible, including S EXTRACTOR setup. 
However, the imaging data used for object detection varies. For 
each single-band fitting run, we have used the same single-band 

^ Note that small systematic sub-pixel offsets between the bands have been 
reported in K12, but we ignore these here. 

^ For reference, the values added for ugrizY J H K -band were 6646, 
2833, 3951, 5610, 25338, 5315, 13620, 59374 and 65015 counts, respec- 
tively. 



Table 1. Summary of the 3 different fitting modes used in this paper 



Mode 


#GALAPAGOS runs 


Method of detection 


Method of titting 


Mode.Sl 


9(1 mn each band) 


individual bands 


single-band 


Mode.S2 


9(1 run each band) 


co-added image 


single-band 


Mode_M 


1 


co-added image 


multi-band 



image for object detection, using an identical SEXTRACTOR setup 
for each (Mode_Sl in Table [T] and throughout this paper). For our 
multi-band fitting run (Mode_M), we used a co-added image that 
was created by simply adding all individual images without further 
normalizatiorLj. Given the background level and noise of the data, 
this creates a slightly biased image towards red bands (e.g. faint, 
u-band-only detections may vanish in the noise from the A'-band 
image), but by adding up the images, we increase the S/N for each 
pixel in the image, allowing for a very deep and detailed object de- 
tection. Most very blue objects, that would not be detected in red 
bands, are bright enough to still be detected in the co-added im- 
age. Again, the same SEXTRACTOR setup is used for this co-added 
image. However, the number of detected objects increases dramat- 
ically compared to all individual bands (see numbers in Table |2]l- 

Using this arrangement, we run the new codes on each single- 
band (ugrizY J H K) individually (Mode_Sl) and one multi-band 
dataset (Mode_M). Bear in mind that, in the multi-band approach, 
only the detection uses a co-added image, the fitting process utilises 
all nine bands individually, though simultaneously. The resulting 
ten catalogues (nine single-band, one multi-band) from these runs 
are then matched using a simple RA/Dec. source correlation and 
all parameters are copied into the one single catalogue that is 
used throughout this section. For consistency checks and to check 
whether pre-existing knowledge about additional neighbouring ob- 
jects gives Mode_M fitting an unfair advantage over Mode_Sl fit- 
ting, we have also run all nine single-band fits using detection on 
co-added images (Mode_S2). The results of this test are discussed 
in §|43]and §|53 

Along with additional information, the resulting catalogue 
contains, for each mode: 

• all SEXTRACTOR output 

• all setup values for GALFIT 

• all fitting values, including values and fitting times 

• all uncertainties for the fitting values 

• file names and folders 

• flags for neighbours, fitting status, constraints hit during the 
fitting process 

• software versions 

Please note that all magnitudes in the catalogue and used 
throughout this paper are total-Sersic magnitudes as directly re- 
turned by GALHT, i.e., integrating the light profiles to infinity. It is 
perhaps more physically meaningful to consider magnitudes based 
on integrating the Sersic profile to some finite radius, e.g. lOro, but 
this is not done in this work. We refer the reader to K12 for further 
discussion of this issue. 

As our sample contains all detected objects instead of objects 
with r < 19.8 for which spectra have been obtained, redshifts are 



^ All GAMA data provided is normalized in the AB system, co-adding the 
images is the correct way of combining the different images in an energy 
driven fashion. 
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only known for the minority of objects in our sample. Requiring 
these would reduce our sample dramatically and restrict the analy- 
sis to a much smaller number of bright objects. Instead, for the ma- 
jority of this paper, we work with apparent magnitudes and sizes, 
giving Tc and other sizes in pixels. This allows us to consider all 
objects measured by GALAPAGOS. Physical values are only used in 
the last sections of this paper where we restrict our consideration to 
those objects with known spectroscopic redshifts, but use a much 
larger area to create a catalogue for a sufficient number of galaxies. 

Before examining the parameters from our GALAPAGOS cata- 
logue, it must first be cleaned in order to select only the objects that 
have been successfully fit by GALFITM. In particular, we wish to 
identify and discard fits with one or more parameters lying on (or 
very close to) a fitting constraint, as described in Section |3] Such 
a fit is unlikely to have found a true minimum in space and is 
indicative of a serious mismatch between the model profile and the 
object in question. This also serves to remove stars from the cata- 
logue, as explained in Section [3] We keep all objects which meet 
the following criteria: 

• minput — 5 < m < mi 

nput + 5, 

• < m < 40, 

• 0.201 < n < 7.99, 

• 0.301 [pix] < rc < 399.0 [pix\. 

• 0.001 < g ^ 1.0, 

• 95 — 5 * magjbest < fwhrajmage and fwhmjmage < 
1 \pix]. These relations were found in the mag-size diagram to 
well separate saturated stars (unsaturated stars) from galaxies and 
is used within GALAPAGOS for this purpose. 

• flag = 2. 

The magnitude input, minput, is the SEXTRACTOR MAG_BEST 
for each object, with offsets where the multi-band detection im- 
age is used. This offset was empirically determined using previous 
results, to adjust (on average) magnitudes measured on the multi- 
band detection image to those for individual bands. The third to 
fifth criteria are slightly more restrictive versions of the fitting con- 
straints used, the next criterion is aimed at separating stars from 
galaxies in the magnitude-size plane, and the last uses a flag that 
is returned by GALAPAGOS. This flag is used to keep track of 
the fitting status of objects and is initially set if the fit has not 
started/tried, 1 if fit has been started ~ e.g. it stays at I if the fit 
failed for some reason- and 2 if GALFIT finished the fit and re- 
turned a result. 

This cleaning is done on a band-by-band basis, i.e., the de- 
cision for each band is entirely independent of the others. For the 
multi-band catalogue, however, we apply these criteria to all bands 
simultaneously, i.e., if values of the fit fail to meet the above crite- 
ria for any band, the entire fit is considered unsuccessful (but not 
the entire polynomial is checked, e.g. interpolated values could in 
places violate these criteria). Hereafter, we refer to the objects in 
this cleaned catalogue as 'successfully fit' or as having a 'valid 
fit result', to distinguish them from objects that were detected by 
SEXTRACTOR, but for which GALFITM 'failed' to find a valid fit 
(flag = 1 and the objects violating the above criteria), and, in the 
case of simulated data, galaxies that were too faint to be detected at 
all. Tables [2] and [3] summarize the numbers of objects and success 
rates. The reader should be advised here to be careful in interpreting 
this success rate as a true success rate (especially when comparing 
to success rates in K12). On the order of half the detected objects in 
SDSS/GAMA imaging are stars, e.g. we would both want and ex- 
pect those to be filtered out by our catalogue cleaning, e.g. lowering 
the success rate to 50%, although every single galaxy could have 



valid fitting results. Whereas KI2 applies techniques to recover fit- 
ting values for objects that initially failed (e.g. by using different 
settings and re-running the fit), such a scheme is not present in our 
software. Although this could be introduced, it creates a less ho- 
mogenous dataset and for the purpose of the single to multi-band 
comparison in the context of MegaMoiph, we avoided this by sim- 
ply ignoring these failed fits in our analysis. 

4.3 Results 

After running GALAPAGOS in both Mode.Sl and Mode_M, clean- 
ing the catalogues and correlating the objects, using the same 
codes, same procedures and with a setup as similar as possible, we 
are now in a position to compare the single and multi-band tech- 
niques. Some example images, fits and fitting residuals for all bands 
are shown in the Appendix. 

Figure[T]shows some fitting results as a function of wavelength 
for three of the objects in our real GAMA sample (top to bottom 
panels). The left column shows recovered magnitudes, as a function 
of wavelength, from both single and multi-band fitting, in compar- 
ison to GAMA photometric data and single-band fits performed by 
K12. The middle column shows the sizes recovered for the same 
galaxies, the right column shows Sersic indices. Please note that 
the X-axis in this figure - and all figures throughout this paper - 
shows linear scaling with wavelength. Although log{X) might be 
physically more meaningful, the scaling parameter in the fitting 
process was chosen to be linear with wavelength in this work and 
this should be resembled in the figures, e.g. a linear polynomial 
would only then appear linear in the figures. Similarly, the slightly 
distorted shape in the plots for size and Sersic results from the log- 
arithmic scaling of the j/-axis. 

One visible effect is that the magnitudes recovered by both fit- 
ting techniques are nearly always brighter than those from aperture 
photometry. This offset is expected, as aperture photometry always 
misses some fraction of the light, whereas the magnitudes from 
GALFIT integrate the profile out to infinit)|fl 

It also becomes clear that even in case of bright galaxies, some 
of the single-band fits (e.g., u- and 2-band in the second exam- 
ple), fail to retuiTi a valid result (with 'valid' being defined in Sec- 
tion l4.2t . For fainter galaxies (e.g. the lowermost example), the suc- 
cess rate for single-band fitting decreases and the scatter increases 
with respect to multi-band and aperture-based results. Bear in mind 
that for magnitudes, we do not constrain the fitting values directly; 
the smoothness of the recovered SED is an indirect result of con- 
straining the profile structural parameters, and not forced by direct 
constraints on the magnitudes themselves. 

For multi-band size measurements, by design, the multi-band 
fitting results lie on smooth curves, which greatly reduces the scat- 
ter in this parameter. Especially in the lowermost example, the 
single-band values vary strongly (and arguably un-physically, with 
a size difference between g- and z-band of nearly a factor of 100). 
Generally, for these relatively bright galaxies (chosen to be GAMA 
spectroscopic targets, with r < 19.8), the sizes from multi-band 
fitting follow the trends of the single-band results, but with a more 
physically realistic smoothness. Even for relatively bright galaxies, 
single-band sizes vary greatly from one band to the next, often by 
a factor of a few. Please keep in mind that the eiTor bars shown in 
all these plots are parameter uncertainties as returned by GALFIT, 

* For a more detailed discussion of this effect o f Sersic profiles, please see 
external literature, e.g. iGraham & Drived fcOOSi) 
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Figure 1. Example fitting results for thi'ee individual galaxies as a function of wavelength. The images, models and fitting residuals of both single-band 
(Mode.Sl) and multi-band fitting (Mode_M) are shown in Figs. lAll to lA3l in the Appendix. We show recovered magnitudes on the left, where we compare 
to both GAMA (SEXTRACTOR MAG_AUTO from catalogue ApMatchedCatv03; green triangles) values and single-band fits by (K12, orange boxes). Black 
asterisks show multi-band results, blue diamonds show results from single-band fits (please note that single-band values have been slightly offset to the right, 
values from K12 slightly to the left, in order to make the error bars more visible). Error bars shown in all panels are parameter uncertainties as returned by 
GALFIT. It is clear that single band fitting introduces larger scatter than multi-band fitting, particularly for fainter objects. In the middle column we show the 
same comparison for galaxy sizes. For comparison, we show the from the Sersic fits of K12. Here, by design, the multi-band fitting returns sizes which 
vary as a polynomial of second order with wavelength, indicated by the black line. In the light column, we show similar figures for Sersic index. Again, by 
design, multi-band fitting returns a smooth dependence over wavelengths. Multi-band fitting generally follows similar trends, but largely reduces the scatter 
and the error bars. In all panels, please note that missing single-band points (both in our own analysis and in the values by K12) indicate that these fits were 
unsuccessful, and so no sensible value is available, e.g. all three of these objects have no single-band fit for the M-band in our own analysis and only the first 
object has a value by K12. Please see the text for more details. 



which have been shown to underestimate the true values (H07) and 
should be interpreted as a lower limit of the true uncertainty. While 
we do not believe these error bars to be realistic, they do allow 
a comparison between single and multi-band fitting. However, we 
would possibly not consider the upturn towards K-band sizes real. 
Especially these galaxies show only individual examples and do not 
represent the population as a whole. An upturn in the K-band sizes 
in the entire population is not found, see e.g. Fig.|5] 



A comparison for Sersic indices is shown in the right column 
of this figure, including a comparison to values of K12 where they 
exist. A comparison to other values from the literature is difficult 
because no such values exist for most of our objects. Generally, a 
trend from lower n in blue bands to higher n in red bands is visible 
for most objects in our sample. We will investigate the recoverabil- 
ity of Sersic indices more in Section|5] where a true value is known 
and an analysis is both easier and more thorough. 
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In the figures of size and Sersic index n versus wavelength 
(middle and right column of Fig.[T}, we not only show the individ- 
ual band sizes for the GALFITM multi-band fit results, but we also 
show the full polynomial function {pi (A; {qi, }), c.f. Section [2.1.1t 
as a black line. An elegant side-effect of fitting these polynomi- 
als, rather than values specific to the wavelength of each band, is 
that it allows easy estimation of sizes (as well as other parame- 
ters) at intermediate wavelengths. For scientific analyses, one often 
wishes to compare restframe parameter values. The polynomial pa- 
rameter functions, inherent to our multi-band approach, provide a 
simple way to determine these, with greater accuracy than could 
be obtained by interpolating between single-band values. However, 
as discussed in Section [S] we cannot take the same approach with 
our magnitudes, because the high-order polynomials used for these 
suffer from Runge's phenomenon. This issue, and ways around it, 
are discussed in Paper I. In this work, we take the conventional 
approach and treat the magnitudes as discrete values, and deter- 
mine restframe values vi a SED template fitting using KCORRECT 
( iBlanton & Roweisll2007h . 

In a similar maimer to the individual examples shown in Fig.[T] 
in Fig.|2]we show trends for magnitudes recovered using both sin- 
gle (both Mode_Sl and Mode_S2) and multi-band fitting for our 
entire sample of real galaxies. Lacking 'true' values, we cannot 
show the offset and scatter of the two methods, instead we consider 
the average SED. All the individual galaxy SEDs were normalized 
to an r-band magnitude of zero before averaging, to minimize the 
scatter due to different galaxy brightnesses. As a comparison, we 
show the average SED for a bright galaxy sample based on aperture 
photometry from the GAMA survey, normalized in the same way. 
This comparison, while not being perfect due to differences in the 
samples shown, gives an indication of the intrinsic scatter in galaxy 
SEDs. 

In the upper panel of Fig.|2] we show the comparison between 
Mode_S I fitting (e.g. only using single-band data for the entire pro- 
cess, including object detection) and multi-band fitting. In the lower 
figure, we show what happens when multi-h&nA detection is used 
for single-hanA- fitting (Mode_S2). 

First, we will discuss the upper figure here. Overall, both sin- 
gle and multi-band fitting show the same trend. Both results show 
slight offsets with respect to the general GAMA SED (as deter- 
mined by aperture photometry, see discussion below), as derived 
from 972 GAMA objects identified in the region. Most of the scat- 
ter in the normalised SEDs is due to intrinsic variation between the 
galaxies. Lacking a 'true' comparison value makes it difficult to 
make more stringent tests, but there are hints that the scatter (for 
the same sample please compare dark-blue and red data points in 
Fig.|2j is slightly reduced in most bands when multi-band fitting is 
used (especially in low S/N bands, u and 2). The normalised SEDs 
for the entire multi-band sample (orange) shows larger scatter and 
offsets even when compared to the full single-band sample (light 
blue). This is a result of the multi-band sample containing fainter 
galaxies than the others. Using real data, it is not clear whether the 
increased scatter compared to the general GAMA SED is a result 
of worse fitting results, or real variations that are not reflected in 
aperture photometry. However, this effect will be examined using 
simulated data without such scatter in Section|5] 

The most dramatic advantage of multi-band fitting becomes 
apparent when comparing the sample sizes of galaxies for which 
parameter values can be derived. Whereas the single-band fits re- 
turn valid results for between 209 (it-band) and 6804 (i-band) 
galaxies, the multi-band fitting returns valid values in all bands for 
15666 objects. However, there are two effects at work here; the 
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Figure 2. The average normalised SED for various methods and samples. 
The upper figure shows Mode.Sl as blue data points, the lower figure 
uses Mode_S2 (increasing sample sizes, but introducing additional scatter), 
please see text for discussion. All input SEDs are normalised to zero magni- 
tude in the r-band. Dark-blue crosses and red asterisks indicate the average 
SED from single and multi-band fits to an overlapping sample, in which 
each galaxy has vaHd fitting results from both techniques. Light-blue dia- 
monds display the same trend for all single-band fits, irrespective of whether 
the multi-band fitting result and especially the other single-band fits in the 
other bands were valid. Conversely, orange triangles show the SED for all 
multi-band fits, in'espective of whether the single-band fit returned a valid 
result. Points other than dark-blue are slightly offset in wavelength for read- 
ability of the figure. The error bars on each point show the rms scatter of the 
normalised SEDs (resistant mean, iteratively clipped at 3 sigma). The num- 
bers above the points for each band refer to the number of valid single-band 
fits for that band (light blue sample). The numbers in the legend give the 
ranges of the other samples. Note that multi-band fits are either valid for all 
bands or for none at all, so the multi-band sample contains the same num- 
ber of objects in each band. For compaiison, the green shaded area shows 
the normalised SED, and its rms scatter, for GAMA aperture photometry of 
972 GAMA objects that were identified in the same survey area, e.g. they 
would be a subset of galaxies in our total sample. As our samples do contain 
fainter objects, not targeted by GAMA, no perfect match is expected. In the 
lower panel, we can see that the improvement of Mode_M versus Mode_S2 
appears very large. However, please see the text for further discussion. 



number of valid fits in each band depends upon both the number 
of objects detected and the fitting success rate. The former is a re- 
sult of the chosen detection image and S EXTRACTOR setup. The 
Mode_Sl results in this figure use single-band detections, while 
the multi-band results are based on detections on a co-added (and 
hence deeper) image. Table |2| gives the number of objects de- 
tected in the imaging for each band, and the number of those ob- 
jects which are successfully fit by single-band GALAPAGOS. It also 
shows the resulting numbers of objects with valid single-band mea- 
surements for every band (ugrizYJHK), or just the six highest 
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S/N bands (griYHK). Finally, Table [2] gives the number of ob- 
jects detected in the co-added multi-band detection image, and the 
number of these with valid multi-band fits (and hence meaningful 
measurements in all ugrizY J H K bands). While the number of 
detections plays an important role (e.g, the number of u-band de- 
tected sources is 20 per cent of that for r-band), the fit success rate 
is significantly higher for multi-band than for any of the individ- 
ual single-band fits, and much greater when one requires complete 
multi-band data. While of course benefiting from detecting more 
objects, this shows that the multi-band approach is more stable and 
thus more often returns valid measurements that are more likely to 
resemble the true parameters of the galaxy. 

The substantial difference in the number of objects returned 
by Mode_S 1 and Mode_M methods is partly due to the initial object 
detection. The multi-band fits are based on a multi-band detection 
image, and hence fits are attempted for many objects which are un- 
detected in some of the single-band images, especially u-band. It is 
possible that single-band fitting may be able to return meaningful 
values for these objects, if it were to be aware of their presence. 
To investigate this, we have repeated the above analysis, but using 
the multi-band detection image even when fitting single-band data 
(Mode_S2). In addition to making the single-band method attempt 
to fit more targets, the additional objects will result in differences in 
deblending, masking and starting parameters. This potentially gives 
single-band fitting a better chance of measuring reliable galaxy pa- 
rameters. 

We will discuss this test in more detail in Section |5] where 
performance comparison is easier. Here we only show an exam- 
ple for magnitude, in the lower panel of Fig. |2] Two effects are 
strikingly evident. Firstly, Mode_S2 fitting now indeed does return 
valid fitting results for many more objects (as apparent from the 
numbers at the top of the figure). The samples are still smaller than 
for multi-band fitting, indicating that GALFIT still fails to return a 
valid fit more often than for multi-band fitting. Secondly, those ob- 
jects that the code does succeed on have a much larger scatter than 
the multi-band results (compare dark blue to red error bars, which 
show the scatter for an identical sample of objects). While the use 
of multi-band detections does greatly increase the sample sizes, the 
scatter in the single-band fitting results has increased dramatically 
as well, especially in the u- and z-bands. This confirms that the im- 
provement in fitting quality and sample size when using multi-band 
fitting does indeed result from the strength of the multi-band fitting 
approach presented in this paper, and not simply due to the advan- 
tage of using a multi-band detection image. Our multi-band fitting 
technique can recover reliable measurements in bands where the 
objects are too faint to be reliably fit (using a single-band method) 
or even detected. 

The relatively poor behaviour of the single-band fits to multi- 
band detections can be attributed to a number of potential causes. 
The simplest explanation is that without the constraints on the size 
and shape of the profile, which naturally come from higher S/N 
bands in multi-band fitting, single-band fits produce more uncer- 
tain results in the low-S/N bands. However, there are two addi- 
tional effects that may be at work. Firstly, given that the fitting po- 
sition of the profiles is only constrained to be within the postage 
stamp, potentially one or more neighbouring (secondary) objects 
may 'wander off, away from their intended targets (which may 
be undetectable in the single band) and on to the (primary) target 
object. Such behaviour would lead to a fainter magnitude being 
returned for the primary, as its flux is distributed between multi- 
ple profiles. Secondly, the opposite is possible. When the primary 
source is much fainter than any of the secondaries, or even invisible 



in the single-band image that is being fit, the primary profile may 
'wander off' to settle on a secondary object. While the stacked pro- 
files would split the secondary flux, the resulting magnitude might 
be brighter than that of the true primary source. We did not further 
investigate which of these effects dominates, as GALAPAGOS does 
not return the fitting values of secondaries and this investigation is 
beyond the scope of this paper. In the case of multi-band fitting, 
these issues becomes much less significant as the profile position is 
effectively constrained using information from all the bands. 

Our Mode_S2 fitting results could hence potentially be im- 
proved by either imposing tighter constraints on the positions, thus 
preventing profiles from 'wandering off their intended target, or 
by using profile information from one fit (e.g., on the r-band) to 
constrain or fix parameters in a subsequent fit on a lower-S/N band. 
However, this is not as natural nor effective a solution as the multi- 
band approach we advocate in this paper. 

The last thing to note from our detection-image test is that 
some of the objects that are missing from our Mode_S 1 fit results, 
but which are recovered in Mode_S2 fits, have relatively bright 
magnitudes. These may be objects with low surface brightness, 
which are undetected in some of the single-band images, particu- 
larly li-band, despite their integrated brightness. Furthermore, even 
when single-band fits are aware of these objects in Mode_S2, these 
fits frequently fail to extract meaningful information, in contrast to 
multi-band fitting. The implication of this is that, without taking 
the multi-band approach, we preferentially lose information about 
galaxies with less-peaky profiles, i.e., disks. 

Figure [3] is similar to Fig. |2] but only considers two 
stringently-defined subsets: 84 galaxies for which all single band 
fits returned a valid result, and 901 galaxies for which we ob- 
tained valid values in the six highest S/N (griY H K) bands. This 
makes the comparison between the codes much easier and cleaner, 
as the same sample of galaxies is compared at all times. Both the 
Mode_Sl and Mode_M fitting results closely agree with one an- 
other and the GAMA aperture photometry. However, while subtle, 
it is also apparent that multi-band fitting slightly reduces the scat- 
ter for most bands. This is shown more clearly using simulations in 
SectionlS^ 

We do not explicitly show a figure similar to Fig. |3] for 
Mode_S2 results. However, we can show from our analysis dis- 
cussed above, that, while again the sample size increases when 
multi-band detection is used, the scatter of the single-band fitting 
results increases dramatically. 

The sample sizes in Fig. [3] highlight an important problem 
with single-band fitting, at least in the simplest case of treating each 
band completely independently, as we do here. As detailed in Ta- 
ble [2] we only obtain meaningful parameter values in all bands for 
87 objects (992 objects in the case that only griYHK values are 
needed). With multi-band fitting, we derive valid parameters in all 
bands for 15666 objects, an increase in sample size of a factor of 
~ 180 (~ 16), and obviously a valuable advantage. As mentioned 
above, this is partly due to the use of a multi-band detection im- 
age (which can only be reliably used with multi-band fitting), and 
partly due to an increase in the fitting success rate, thanks to the 
stability achieved by constraining the profile shape across multiple 
bands. The advantage of the multi-band approach becomes espe- 
cially evident in science cases where parameter values are required 
in several bands, which, in case of single-band fitting, decreases the 
sample size dramatically. While multi-band fits derive valid values 
for ~ 54 per cent of detected objects (~ 71 per cent of galaxies 
considering only a sample of bright GAMA galaxies), a combined 
single-band approach only recovers valid parameters for < 0.4 per 
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Table 2. The number of objects in the real dataset that are detected and successfully fitted, and the corresponding success rate, for each single-band image 
(only Mode.Sl), for the Mode.Sl results combined, and for the multi-band (Mode_M) technique. 'Combined' in this context does not mean detection on the 
co-added image (Mode_S2), but combination of all the single-band detections and fit results to create a sample of all detections and successful fits when using 
ugrizY J H K -hanAfi or griY H K -hands,, respectively. The reader should take care in interpreting this success rate as a true success rate. On the order of 50 
per cent detected objects in SDSS/GAMA imaging are stars, that are thrown out during the cleaning of the catalogue. This lowers the success rate dramatically, 
even if all galaxies had valid fitting results. An interpretation of these rates in an absolute sense should take this effect into account, but is not given in this 
analysis. The fraction of stars in the detections is possibly higher in the shallower bands due to their higher detectability, which would partly explain the lower 
success rates in these band. As a reliable star/galaxy classification is difficult in order to correct these numbers, we instead show the same values on the right 
side, but restrict the sample to targeted GAMA galaxies, e.g. no stars and only galaxies with r < 19.8. Success rates are indeed much higher when only this 
bright galaxy sample is considered. Comparing the numbers at the bottom of the table, we can successfully derive all fitting values for 885 GAMA objects out 
of the 1251 GAMA objects that were detected when using multi-band fitting. Requiring all ugrizYJHK fitting results from single-band fitting reduces this 
sample to 75 objects. 



Band 


all objects 

#objects fit #objects detected Success rate 


GAMA objects 
# fit # detected Success rate 


u 


209 


3773 


5.5% 


108 


281 


38.4% 


g 


3276 


13547 


24.2% 


827 


1229 


67.3% 


r 


5802 


19169 


30.3% 


896 


1237 


72.4% 


i 


6804 


22059 


30.8% 


907 


1235 


73.4% 


z 


2117 


12496 


16.9% 


725 


1178 


61.5% 


Y 


3912 


16200 


24.1% 


859 


1221 


70.4% 


J 


2530 


13289 


19.0% 


686 


1150 


59.7% 


H 


5432 


17750 


30.6% 


885 


1212 


73.0% 


K 


5190 


16117 


32.2% 


770 


1171 


65.8% 


combined single-band (ugrizYJHK) 


87 


24995 


0.35% 


75 


1251 


6.0% 


combined single-band (griY H K) 
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cent of these objects (6.0% of bright galaxies), and ~ 4 (^^^ 44) 
per cent when only values in the higher S/N bands {griY H K) are 
needed. Such an science case is discussed in Section|7] 

The normalised SEDs, given in Figs.|2]and|3]allow us to eval- 
uate the band-to-band variation of our measurements, and compare 
these to aperture photometry. However, the scatter in these SEDs 
are dominated by intrinsic variations in the galaxy population. We 
therefore perform a further comparison of our GALFIT magnitudes 
directly against GAMA aperture photometry, a test that would be 
carried out by most authors as performance tests of profile fitting 
codes. The results of this test are shown in Fig. Ill 

This comparison in Fig. |4] is only possible for a subset of 
galaxies in each band, those which have valid fitting results for 
the method in question. As previously mentioned, some offset is 
expected as aperture photometry will always miss some fraction 
of the light from the outskirts of the galaxies. On the other hand 
integrating galaxy light profiles out to infinity can add fictional 
flux particularly for high-n systems, making the definition of a true 
value difficult. The GALFIT magnitudes are indeed systematically 
brighter than aperture photometry by 0.1-0.2 mag. The scatter be- 
tween GALFIT and aperture magnitudes is also 0.1-0.2 mag, de- 
pending on typical S/N of the band, suggesting that much of this 
scatter could result from variations in the flux missed by aperture 
photometry for galaxies of different profile shapes. The single and 
multi-band fitting results appear quite similar, although the multi- 
band magnitudes generally display slightly greater offsets. How- 
ever, from this figure, it is difficult to determine which magnitudes 
best represent the total light and colours of the target galaxies. We 
explore this in a more satisfactory manner in Section |5] with the 
use of simulations. 

Magnitudes are integrated quantities, and are therefore gener- 



^ This figure looks very similar if we use magnitudes from the individual 
SEXTRACTOR rans peribmied for object detection in GALAPAGOS. 



ally the easiest galaxy properties to extract from imaging. Another 
important galaxy property is size. Size is a natural, direct product 
of Sersic profile fitting. It may be defined in many ways, but a con- 
venient and commonly used value is the half-light radius. This is 
the size measure returned by GALFIT. 

We follow the same approach as for magnitudes (c.f. Figs.|2] 
and[3j. Figures |5] and |6] show the trend in galaxy half-light radii, 
normalised to the r-band, a s a func tion of wavelength (similar 
trends have been reported by iKelvine t al. 2012, see Fig. |6l. Error 
bars indicate the scatter in these normalised sizes. Figure [5] com- 
pares the Mode.Sl and Mode_M results using samples of valid ob- 
jects for each technique, and an overlap sample, containing only 
objects with valid results in both techniques (dark blue and red dat- 
apoints). Figure |6] limits the comparison sample further to those 
for which single-band fits for all bands, or the six highest S/N 
(griY H K) bands, returned a valid result. In contrast to magni- 
tudes, for which the scatter is only slightly improved, the scatter 
on the recovered sizes are hugely improved by multi-band fitting. 
These figures also display an average trend, recovered by both tech- 
niques, such that galaxies appear smaller in redder bands. This can 
be explained by two processes. Firstly, blue bands galaxy light is 
typically dominated by galaxy disks, whereas in red bands bulges, 
which are usually smaller than disks, contribute a greater fraction of 
the light (see magenta and orange lines in Fig.|6l. Secondly, a com- 
bination of dust and s tellar populations in pure disk galaxies can 
create a similar effect dMollenhoff et al.ll2006t|Pastrav et alj|2012l; 
Pastrav et al., in prep., see dashed line in Fig. |6]for the effects of 
dust on the measured half-light radius). 

The reduction in scatter we see in Figs. |5]and|6]is at least 
partly by design, as the half-light radii are constrained to lie on a 
second order polynomial as a function of wavelength. Such a re- 
duction in scatter would therefore be seen even if the real sizes of 
galaxies did vary wildly with wavelength. However, from inspec- 
tion of the single-band fits, i.e. figures like Fig. [T] for many more 
examples, we see that sizes do show a smooth, and nearly linear. 
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Figure 3. Similar to Fig.|2] tliese figures sliow the average normalised-SED 
and the scatter around that average. Both figures use Mode.Sl data, but now 
for a small sample of 84 galaxies for which all single-band and multi-band 
fits return a valid results (top panel) and a sample of 901 galaxies for which 
the griYHK bands return valid results (bottom panel). All points within 
each figure therefore incorporate the same objects (unless numbers indicate 
differently in the lower panel). The shaded green areas in each figure show 
normalised-SEDs based on GAMA aperture photometry, using the galaxies 
in our samples for which these values exist. This sample hence forms a very 
similar sample to the ones used in the plotting symbols, which explains the 
much better agreement between points and shaded green area compared to 
FigE 



trend with wavelength. Consideration of physically realistic stel- 
lar populations also supports smooth variations in half-light radius 
with wavelength. Such was our justification for using second or- 
der polynomials in the first place. As we lack an independent esti- 
mate for the 'true' sizes of the objects, we can not investigate this 
further with real data. However, in the following section we use 
simulations to demonstrate that the multi-band approach results in 
significantly improved size measurements, suggesting that the re- 
duced scatter in Figs.|5]and|6]corresponds more closely to the true 
variation in relative size with wavelength. 

For comparison, we illustrate the trends described in K12 for 
both spheroid-dominated and disk-dominated galaxies in Fig. [6] 
Given the clean sample in this figure, which mostly consists of 
bright galaxies, that would also be targeted by K12, we can con- 
firm to find very similar trends of with wavelength. 

Determining galaxy size may be thought of as finding a radius 
for which the integrated flux within and without corresponds to a 
given ratio. The behaviour of the profile within each of those two 
regions does not matter, and there is no ambiguity in the definition 
for any arbitrary monotonic profile. Determining Sersic index re- 
quires further information, particularly regarding the behaviour of 
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Figure 4. The difference between magnitudes obtained from GALFIT (blue 
crosses: Mode_Sl, red asterisks: Mode_M) and aperture photometry mea- 
sured by the GAMA survey, using SEXTRACTOR. We show the figure for a 
sample of galaxies with r < 19.8 for which, for each filter band, all of these 
measurements exist. We recover fainter magnitudes below the green line, 
and brighter magnitudes above. At first glance, the multi-band fitting results 
seem slightly worse, returning significantly brighter magnitudes. We will 
demonstrate, however, in Section |5] specifically Fig. 1121 that SEXTRAC- 
TOR magnitudes are bad proxies for this kind of test and that the multi-band 
fits are better at recovering true total magnitudes. 
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Figure 5. A comparison of the trend in measured half-light radius versus 
wavelength for single and multi-band fitting. Symbols and colours indicate 
the same samples as in Fig. |2] The half-light radii for each galaxy are nor- 
malized using the r-band, and the average variation in relative size versus 
wavelength is shown by the points. We recover larger half-light radii above 
the line, and smaller below. The rms scatter around this trend is indicated by 
the eiTor-bars. Both single and multi-band approaches reproduce no signifi- 
cant average trends, but clearly the multi-band technique reduces the scatter 
substantially. 



the profile at its peak and in its tail. Higher Sersic indices imply a 
greater shift of flux from around the half-light radius and into both 
the centre and outskirts of the profile. However, these pieces of in- 
formation need not be consistent for an arbitrary monotonic profile. 
The Sersic index is therefore more difficult to measure accurately 
and consistently. 

In Fig. [7] and [8] we present a comparison of Mode_Sl single 
band and Mode_M multi-band fitting results for Sersic index along 
the same lines as for half-light radii, again, slightly arbitrarily, nor- 
malizing all objects to the r-band Sersic index. Again, a general 
trend in n with wavelength, similar to the ones reported in K12, is 
recovered by both single and multi-band fitting, although, given the 
scatter in single-band, it is much more visible in multi-band results. 
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Figure 6. Similar to Fig. [5] these figures compare trends of half-Hght radius 
versus wavelength, but now for (top panel) a small sample of 87 galaxies 
for which all Mode.S 1 and Mode.S 1 fits return a valid results and (bottom 
panel) a sample of 901 galaxies for which the griYHK bands return valid 
results in Mode.S 1 . All the points within each figure therefore correspond 
to the same objects. The improvement in scatter when using multi-band 
fitting is especially evident in this figure. For comparison we show the rela- 
tions described in K12 for spheroid-dominated galaxies (orange) and disk- 
dominated galaxies (light blue) which describe the same trends found in our 
data. The black dashed line shows the effect of dust content in the galaxy 
on the size measurements (while ignoring intrinsic stellar population gradi- 
ents) from Pastrav et al (in prep), based on the radiation transfer model of 
|Popescuetai]i201lh . 



This trend confirms our argument that in blue bands one fits a light 
profile dominated by a galaxy disk (low n) whereas in red bands, 
the bulge dominates the light profile (high n). The same effects re- 
garding the scatter of the fitting values are seen, multi-band fitting 
hugely reduces the scatter in n and increases fitting quality. As ex- 
pected, the scatter in both single-and multi-band fitting is larger for 
n than for r^, strengthening the argument that Sersic indices are 
indeed the hardest values to measure. 

We can compare our data to the trends described in K12, 
which we show as lines in Fig. [8] We generally find a qualita- 
tively similar trend, but find it to be stronger than discovered in 
K12. The reason for this is that the sample in K12 has been split 
up into disk-dominated and spheroid-dominated samples using a 
cut in the (U-V)-nfc plane. This split and the normalization to the 
r-band Sersic index in Fig.[8]suppress the trend found in our work 
and lead to shallower trends. This effect can be confirmed in Fig.[8] 
As splitting up the sample into disk and spheroid dominated sam- 
ple is beyond the scope of this paper, this effect is expected in this 
comparison. 

In a similar manner to half-light radius, the reduction in scatter 
for Sersic index is partly by design, as we constrain the values to 
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Figure 7. A comparison of the trend in Sersic index versus wavelength 
for Mode.Sl and Mode_M fitting. This figure includes all galaxies, e.g. no 
'morphology' selection has been done. Symbols and colours indicate the 
same samples as in Fig. [2] The Sersic indices for each galaxy are normalized 
using the r-band, and the average variation in n versus wavelength is shown 
by the points. We recover larger Sersic indices above the fine, and smaller 
below. The rms scatter around the observed trend is indicated by the en'or- 
bars. Single and multi-band approaches reproduce similar average trends, 
but the multi-band technique reduces the scatter substantially. 
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Figure 8. Similar to Fig.|3] these figures compare trends of Sersic index ver- 
sus wavelength, but for a small sample of 87 galaxies for which all single- 
band and multi-band fits return a valid results (top panel) and a sample of 
901 galaxies for which the griYHK bands return valid results (bottom 
panel). All the points within each figure therefore coiTcspond to the same 
objects. The improvement in scatter when using multi-band fitting is very 
evident in this figure. 
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lie on a second order polynomial as a function of wavelength with 
the caveats discussed above. Again, from inspection of the single- 
band result, we decided that fitting a second order polynomial to the 
Sersic index as a function of wavelength generally provides suffi- 
cient flexibility to fit most real objects. Again, we refer to 15. 21 for 
a more detailed and normalization-independent evaluation of the 
fitting accuracy of the Sersic parameter. 



5 APPLICATION TO SIMULATED IMAGING 

Unlike observed data, which has the advantage that the galaxies in 
the data truly are real and include internal galaxy structures, such as 
spiral arms, star-forming regions and asymmetries, simulated data 
has the advantage that one knows the input galaxy parameters. This 
allows for a direct comparison of recovered magnitudes, half-light 
radii and Sersic indices to their true values, instead of being limited 
to comparing to alternative estimators, such as aperture photome- 
try or Petrosian radii, or similar Sersic profile measurements from 
other studies. The results are therefore much cleaner and free of 
the intrinsic scatter which dominates most figures in the previous 
section. Additionally, one can examine in more detail the ability of 
multi-band fitting to recover the variation of galaxy sizes and Sersic 
indices with wavelength. 

On the other hand, simulations have the disadvantage that the 
objects now all possess the smooth, symmetric profiles for which 
we know that Sersic profiles (with some polynomial wavelength 
dependence of their parameters, and convolved by a known PSF) 
will precisely fit the data. Our simulation can only test the fitting 
routines under the assumption that galaxies are intrinsically repre- 
sentable by perfect Sersic profiles. We know that for some galaxies, 
and for all galaxies observed in sufficient detail, this is not true. In 
this sense, simulations are an idealised case and, while affording 
a good comparison of single-band to multi-band fitting, they can 
only give a lower limit on the error bars for real data. 



5.1 Creating the simulations 

To create our simulated data, we follow the methodology of H07, 
making adaptations to the scripts as necessary to enable the simu- 
lation of a multi-wavelength dataset. 

As in H07, we started from a catalogue of results from fitting 
profiles to real data, in order to get the distribution of galaxy pa- 
rameters in the simulated images as close to that of a real sample 
as possible. We decided to start from the sample used in Section 
|4]and 'cleaned' as described above. We first divide the r-band m- 
log(rc) plane into bins of width 0.4 in m and 0.1 in log(re). Fol- 
lowing H07, to create the parameters for each simulated galaxy we 
select an r-band m-rc bin, and use the centre of this bin as our 
starting value for r-band magnitude. We then pick a real galaxy at 
random from that bin, from which we obtain realistic values for r- 
band effective radius, Sersic index and axis ratio, as well as the 
wavelength-variation for each parameter, including colours. The 
centre position and position angle were chosen randomly for each 
object. 

To produce a simulated dataset with the same magnitude and 
effective radius distribution as the real data, we selected r-band m- 
ro bins with a probability weighted by the distribution of real galax- 
ies in this plane. However, in order to test the code at the detection 
limits we extend the faint side of the magnitude distribution. This 
was achieved by finding the peak in the number density of galaxies 



versus magnitude, for each bin, and extending this peak level 
two magnitudes fainter (as was done in H07). 

As in H07, we slightly varied the galaxy parameters produced 
by the above procedure, in order to create a smooth distribution, in- 
stead of simply re-simulating the same few galaxies over and over. 
These applied variations are as follows: 

Position: The x and y centre coordinates in the simulated field 
are randomly chosen (although of course neighbouring tiles are 
consistent in their overlap regions), but are identical for all bands 
(no statistical or systematic offsets are applied). 

Magnitude: Starting from the initially chosen r-band magni- 
tude, we first apply an offset, drawn from a Gaussian distribution 
with standard deviation of 0.4 mag. To define the magnitudes in 
the other bands, colours (as differences of the other bands to r- 
band) are directly taken from a real galaxy within the corresponding 
m-Tc bin and added to the r-band magnitude. The magnitudes in 
the individual bands are finally modified by adding Gaussian noise 
with a 0.1 standard deviation. The magnitudes were allowed to vary 
freely as function of wavelength in the fitting process applied to the 
real galaxies, so there are no smoothness constraints on magnitude 
in the simulated dataset. 

Half-light radius: The fits to the real data were caiiied out us- 
ing a polynomial of second order to describe the wavelength de- 
pendence of re. We start with the function taken from a real galaxy 
in the selected m-re bin. We then add an additional linear slope, 
drawn from a Gaussian distribution with standard deviation corre- 
sponding to a 10 per cent change in over the full wavelength 
range. No noise is added to the individual bands, so that simulated 
values follow a second order polynomial exactly. 

Sersic indices: Similar to r^, in the fits to the real data the 
wavelength dependence of n was described by a quadratic. For 
most of the simulated galaxies, starting values are taken from the 
same real galaxy as before, but in this case we simply modify the 
values in all bands by a factor drawn from a Gaussian with unit 
mean and standard deviation of 0.3. This offsets the values, but 
preserves the linear and quadratic coefficients. The simulated val- 
ues therefore follow a second order polynomial exactly. However, 
for 5 per cent of objects we randomly chose n uniformly from the 
range 0.2-8, with a linear slope as a function of wavelength. This 
was done in order to cover areas of parameter space that are not 
covered by the cleaned catalogue of real galaxies, either because 
no such galaxies exist or such galaxies defy reliable detection or 
modelling. As these areas are of particular interest to us, we simu- 
late some galaxies to cover them. 

Axis ratio: This is simply taken from the chosen real galaxy 
and offset by a Gaussian random variable with standard deviation 
0.1 (subject to the allowed 0-1 range of this variable). The same 
value is adopted by all bands, and no additional noise is added. 

Position angle: This is randomly chosen, with the same value 
applied to all bands, with no further noise added. 

Following H07, we added Poisson noise to the simulated im- 
ages and then add them into an image that is made up of empty sky 
patches from the real GAMA imaging. This means that sky noise is 
present in the simulated images with the same properties as in the 
real data. 

The result of this procedure is a realistic-looking set of images 
with thousands of galaxies that mostly show the same parameter 
distributions as the real galaxies, and for which we know the true 
parameter values. It should be noted that these images do not con- 
tain any stars, due to the cleaning of the catalogue prior to parame- 
ter selection. A fraction of real galaxy images are contaminated by 
stars, which must be masked or deblended by GALAPAGOS, poten- 
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tially leading to additional uncertainties which are neglected by our 
simulations. However, blended galaxies are naturally included by 
our simulation method, and so the effects of blending are included 
in our results. 



5.2 Results 

In this section, we examine the single and multi-band fitting results 
by comparing the recovered parameter values to the true, simulated 
values directly. Where possible, we recreate figures similar to those 
in Section |4l but we are able to examine issues in more detail and 
add a set of figures which could not be created using real data. 

When running the codes we use exactly the same setup and 
procedures as were used for real data in the previous section. In par- 
ticular, we use the same versions of our codes (GALFITM version 
0.1.2.1, GALAPAGOS version 2.0.2). We run the same three scenar- 
ios, i.e., Mode_Sl, Mode_S2 and Mode_M, and carry out the same 
procedure to combine the individual catalogues (9 single-band, 1 
multi-band) in order to create one single catalogue which is used 
throughout this section. This catalogue is additionally correlated 
to the simulated catalogue using RA and Dec, so that true values 
for all objects are known (disregarding the potential for rare mis- 
identifications). 

We show some fits of individual simulated galaxies in Fig. [9] 
As for real galaxies, it is apparent that multi-band fitting returns 
values for more bands, making studies using multiple sets of val- 
ues (e.g., colours) possible for a larger sample of galaxies. It also 
reduces the scatter on all the values that are measured. Whereas the 
single-band magnitudes (and even sizes in the brighter galaxies) 
generally return sensible values, the and Sersic indices recov- 
ered are typically much noisier in single-band fitting compared to 
multi-band fitting. Multi-band fitting generally recovers the shape 
of the wavelength-dependence well, but this is of course a function 
of magnitude, where fainter galaxies are generally less well recov- 
ered. 

In addition to returning more accurate values at the specific 
wavelengths of the input images, multi-band fitting also provides 
the entire polynomial over the wavelength region. Estimations of 
values at intermediate wavelengths, e.g. restframe values, therefore 
become trivial. As long as they are within the wavelength range 
covered by the input images, they can simply be obtained from the 
polynomial. As mentioned in Section [3] this interpolation is not 
sensible for magnitudes, for which we use a higher-order polyno- 
mial, due to Runge's phenomenon. In any case, we have a better 
understanding of galaxy SEDs, and so it is wise to use this knowl- 
edge, e.g., via template spectra, to estimate accurate magnitude k- 
corrections. 

After visually inspecting a large fraction of individual fits, we 
found that multi-band fitting was able to successfully outperform 
single-band fitting in all cases when the multi-band fit was not re- 
jected due to proximity to a fitting constraint. The opposite case - 
objects where multi-band fails but single-band fits return results for 
at least some bands - does occur, but is less common. Numbers to 
illustrate this are included in Fig. [13] As expected we also find that 
for very faint galaxies (r-band magnitude ~ 23, close to the SDSS 
detection limit), the wavelength dependence of the sizes and Sersic 
indices do not match that simulated, and the magnitudes are very 
noisy. On most of these objects single-band fitting is not able to re- 
cover many, if any, values. However, parameter estimates from the 
Mode_M fits, though noisy, were generally still in agreement with 
the simulated values, given the big error bars. 

In Fig. [To] we illustrate the performance of single- and multi- 
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Figure 10. Magnitude recoverability for simulated galaxies (Mode_M and 
Mode.Sl in the top panel, Mode_M and Mode_S2 in the bottom panel). 
Note especially the change of quantity on the y-axis compared to Fig. |2] 
for simulated galaxies we are able to demonstrate the difference between 
simulated and measured values. This means that the measurement uncer- 
tainties are no longer convolved with the intrinsic scatter on the colours, 
permitting a more stringent comparison between techniques . For example, 
while on real galaxies the scatter in M-band is around 0.5 mag, it is only 
about 0.2 mag here. Multi-band fitting produces a smaller scatter in com- 
parison to single-band results. Systematics are very small, although both 
techniques do display a bias to recover brighter magnitudes by < 0.05 mag. 
The sample size is enlarged by the multi-band approach from 714 (u-band) 
or 12602 (ii'-band) in single-band fitting to 23632 objects in multi-band fit- 
ting, nearly a factor of 2 more objects even for the best single bands. These 
numbers are compai'ed in more detail in table |3). In the bottom panel, the 
bias introduced by combining multi-band detection with single-band fitting 
(Mode_S2) becomes very apparent. This is simply a result of attempting to 
model objects that are below the noise limit. 



band techniques for recovering the magnitudes of our simulated 
galaxies. This is an equivalent figure to Fig. (2] but compares di- 
rectly to the simulated parameter values, rather than arbitrarily nor- 
malizing the SED to the r-band. The intrinsic scatter in galaxy 
colours, which obscured improvements in scatter in Fig. |2] is not 
present in this figure. The offset and scatter in this figure represent 
the true bias and uncertainty in our ability to recover object pa- 
rameters. As in Fig. |2] one can see that the scatter is significantly 
reduced in most bands when using multi-band fitting. Multi-band 
fitting also somewhat reduces the small offset that is apparent in 
low-S/N bands (mainly u- and 2;-band). 

The scatter and bias for the complete multi-band sample 
(shown in yellow) are significantly larger. This is largely due to 
the inclusion of large numbers of faint objects in this sample and 
will be examined in detail later (in Fig. 1 13). 

The lower panel in Fig. [TO] shows the same figure as the up- 
per panel, but for Mode_S2 results. Similar to what was seen using 
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Figure 9. Fits for individual simulated galaxies, similar to Fig. [T] in case of real galaxies. The images of these objects, together with models and fitting 
residuals from both Mode.Sl and Mode_M fitting, are shown in Figs. lA4l to lA7l in the Appendix. The left column shows the magnitude difference between 
fitted and simulated values. The green Une/values indicates an ideal fit (no offset between fit and simulated values), blue diamonds show Mode.Sl fitting 
results, black asterisks show Mode_M results. The middle column shows both simulated and recovered half-light radii in all bands. The black line shows the 
entire polynomial that is fit to the data for Mode_M; in contrast to single-band fitting, where only some of the individual values exist, we derive the entire 
function in the case of multi-band fitting. The right column shows the simulated and fitting values for Sersic indices. The r-band magnitude of the objects is 
indicated in the leftmost figure. See further discussion in the text. 



© 2012 RAS, MNRAS OOO.rTPT] 



20 Hciufiler et al. 



Table 4. Mean and scatter of the magnitude offsets for Mode.Sl and 
Mode_M for a sample of objects with complete Mode.Sl measurements, 
as shown in the top panel of Figure [TT] 



Band 


offset ±a Mode.Sl 


offset ±(j Mode_M 


u 


-0.06±0.17 


0.00±0.09 


g 


0.02±0.04 


0.01 ±0.04 


r 


0.01±0.03 


0.01±0.03 


i 


0.00±0.04 


-0.01 ±0.04 


z 


-0.05±0.12 


-0.02±0.12 


Y 


-0.03±0.11 


-0.03±0.07 


J 


-0.03±0.14 


-0.03±0.09 


H 


-0.02±0.11 


-0.05±0.11 


K 


-0.03±0.13 


-0.04±0.17 



real galaxies, the improvement from Mode_S2 to Mode_M fitting 
is more dramatic thian from Mode_S 1 to Mode_M. Our arguments 
for this are similar to the ones presented in Section |43] The results 
of Mode_S2 fitting are worse largely because the effect of neigh- 
bouring galaxies. Either the primary object wanders off to fit a sec- 
ondary (in which case the fit would likely appear to be too bright), 
or a secondary can centre on the primary (in which case the fit 
would appear to be too faint). Multi-band fitting is largely immune 
to such issues, as the higher S/N images naturally constrain the po- 
sition of the target in the lower S/N bands. 

Table [3] presents the fitting success rates for the simulated 
galaxy sample, in a similar manner to those shown in Table |2] for 
the real data. The fitting success rate here is much higher than in 
case of real galaxies. There are two distinct effects responsible for 
this. Firstly, the simulated images do not contain stars, due to the 
cleaning of the fitting catalogue prior to parameter selection. As 
stars typically result in 'failed' fits, by design, their absence from 
the simulations naturally means that the resulting success rate is 
boosted. The remaining improvement is largely due to the fact that 
all the simulated galaxies actually resemble smooth Sersic profiles, 
which makes the fitting process much more stable, even in single- 
band fitting and in low S/N bands. 

The multi-band technique is not the one with the highest suc- 
cess rate on simulated data. Due to a more sensitive detection in 
Mode_M compared to Mode_Sl, the sample includes many more 
faint galaxies which are harder to fit and decrease the success rate. 
The overall number of objects for which parameters can be de- 
rived is still higher. When comparing the overlap of the samples 
of single-band with each other and the sample in multi-band fit- 
ting, the increase is still substantial, effectively providing fits for 
seven times more galaxies in a sample that is only twice the size to 
start with. This emphasises how dramatically the size of scientific 
samples can be increased by using multi-band fitting. 

To emphasise the effect on sample sizes, and ensure a consis- 
tent comparison. Fig. [TT] again shows the recoverability of magni- 
tude by both- single and multi-band fitting, but for matched sub- 
samples for which all, or most, Mode_S 1 results return valid fitting 
values. In both cases, the reductions in both scatter and offsets are 
clearly visible for most bands. Low S/N bands particularly bene- 
fit from the multi-band approach. The values from the top panel in 
Fig. [TT] are also given in Table [4] Biases are somewhat reduced at 
most wavelengths, while scatter is reduced significantly in the low 
S/N bands and slightly in the other bands. 

In Fig. [4] we compared magnitudes derived by GALFIT with 
those obtained using aperture photometry, finding a consistent off- 
set of ~ 0.1 mag, which appeared greater when using the multi- 
band technique. In Fig.[T2]we follow up on this by comparing mag- 
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Figure 11. Magnitude recoverability in simulated data for galaxies with 
complete sets of single-band fits. These figures are the same as the top 
panel of Fig llOl but show only a common sample of 297 galaxies for which 
all Mode.Sl fits return a result (top) or 2351 galaxies that have fits in 
griYHK (bottom). For these common samples, multi-band fitting clearly 
increases the accuracy and precision for magnitude measurements. For ref- 
erence, the values shown in the top panel are given in Tablei] 
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Figure 12. Equivalent to Fig. |4]we show here the deviations of recovered 
magnitudes for both Mode_M (red) and Mode.Sl fits (blue) compared to 
SEXTR ACTOR 'Mag-Best'. For comparison, we show the deviations of the 
SEXTRACTOR magnitudes from the simulated values as black data points. 
Clearly visible, SEXTRACTOR magnitudes produce a much larger offset to 
simulated magnitudes and are not a good proxy for 'real' magnitude. Either 
multi- or single-band fitting recover the simulated magnitudes better This 
puts Fig.|4]and the offsets of Mode.Sl and Mode31 values in this figure in 
perspective. 
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Table 3. Object numbers in the simulated dataset and fraction with successful fits for Mode.S 1 and ModeJM. Please note 
that the success rate in this simulated sample resembles an actual success rate and does not include stars that are purposefully 
excluded from the final analysis, but which suppress the apparent success rate in Table|2] In that respect, the success rates in 
this table are more comparable to those presented for the GAMA galaxy sample in Table|2] 



Band 


#objects fitted 


#objects detected 


Success rate 


u 


714 


1213 


58.9% 


rr 

5 


6022 


8839 


68.1% 


r 


10362 


13739 


75.4% 


i 


12326 


16490 


74.7% 


z 


4139 


6196 


66.8% 


Y 


9627 


13123 


73.4% 


J 


7965 


11068 


72.0% 


H 


12593 


16992 


74.1% 


K 


12602 


17755 


71.0% 


combined single band (ugrizYJHK) 


305 


26308 


1.2% 


combined single band (griY H K) 


2560 


25631 


10.0% 


mwl 


23632 


31798 


74.3% 



nitudes from Mode.Sl and Mode.M fits, as well as SEXTRAC- 
TOR 'Mag-Best' aperture magnitudes, to their simulated values. 
The GALFIT results are the same as those shown in the top panel 
of Fig. [To] It is obvious that the offset in Fig.|4]is a result of SEx- 
TRACTOR underestimating the total flux. Both single and multi- 
band fitting provide better estimates of total magnitude. Fig. |4] 
therefore does not show which method recovers the better values. 
However, as Fig. [4] together with Fig.[lO] serves as a simple illus- 
tration of the issues in comparing aperture and model photometry, 
we retain it. Note that our simulated magnitudes are based on Sersic 
models integrated to infinity, as are the magnitudes recovered by 
GALFIT. Aperture magnitudes necessarily do not measure the flux 
in the distant tails of the profile, and it is arguable whether this flux 
actually exists, as much of it is at surface brightnesses fainter than 
can possibly be measured. This is discussed in much more detail 
by K12, who advocate using magnitudes obtained from truncated 
Sersic profiles for all science analyses. 

A more complete presentation of the results from fit- 
ting simulated data is presented in Fig. [13] For each band 
(it, g, r, i, z, Y, J, H, K from top to bottom) we show the magni- 
tude offset from the fitting values, but instead of putting everything 
into one bin as in the previous figures, we show them as a function 
of simulated (true) galaxy magnitudj^ The left column shows re- 
sults from Mode.S 1 fits with the mean offset and rms scatter over- 
laid in light blue (dotted lines indicate the ±lcr scatter); the right 
column shows the same for Mode_M fits, with the mean offset and 
rms scatter for the whole object sample in light blue. Black num- 
bers in the top right corner of each panel show the number of ob- 
jects appearing in the figure, i.e., simulated objects which have been 
detected and successfully fit by the method in question. Addition- 
ally, in both columns, we overplot the sample of objects for which 
both Mode_S 1 and Mode_M fits return a valid result in yellow. Yel- 
low numbers in the top right corner of the left-column panels give 
the number of objects in these samples. One can see that not all ob- 
jects with single-band fits have successful multi-band fits, but the 
fraction is much higher than the other way around. We show the 
mean trends and rms scatter for this overlap sample for Mode_Sl 
fits in blue in the left column and, to aid comparison, for Mode_M 
fits in red in both columns. 



^" We note that when using observed surface brightness on the x-axis, the 
figures look qualitatively very similar 



For an ideal fitting code, the mean and scatter lines should be 
close to zero at all input magnitudes, indicating that no systematic 
or statistical biases exist in the fitting routines. As one can see, both 
single-band and multi-band fitting values are not significantly bi- 
ased for most galaxies, but the multi-wavelength fitting reduces the 
scatter somewhat, particularly in lower S/N bands. The improve- 
ment is not significant for the overlap sample of galaxies in some 
higher S/N bands, but even in these cases the multi-band approach 
significantly improves the completeness of the sample. 

At this point, we would like to remind the reader that be- 
fore creating these figures we cleaned our catalogue of objects that 
ended up close to fitting constraints. As mentioned above, the total 
number of successfully fit objects is indicated by the black num- 
ber in the comer of each panel. As one can see, the numbers of 
objects for which we can get valid results is much larger in the 
case of multi-band fitting, confirming our findings from Section|4] 
Whereas in r- and z-band the sample size only doubles, we can ob- 
tain fitting parameters for 30 times more galaxies in the ii-band. 
In most bands, we can also recover values for fainter galaxies. As 
these are harder to fit and result in larger error bars, they al so lead 
to a large 'scatter' in Figs. |10|toll2l(orange symbols). This is ex- 
pected and does not imply that Mode_M fits are less accurate than 
Mode_S 1 fits, as could be wrongly inteipreted from the error bars 
in those figures. 

For reference, we show, as a vertical purple line in the r-band 
panels in Fig.[T3] the m,. = 19.8 magnitude limit of the spectro- 
scopic survey carried out as part of the GAMA project. For galaxies 
brighter than this limit, both single-band fitting and multi-band fit- 
ting do a similar job of returning good fitting values. However, the 
accuracy and precision of the recovered parameters are clearly im- 
proved by the multi-band method, particularly for lower S/N bands 
and, as we shall see, size and Sersic index. One distinct advantage 
of the multi-band technique is its ability to obtain valid results for 
fainter objects. However, if one is only interested in bright galaxies, 
e.g. for which spectroscopy has been obtained, then the difference 
in terms of sample size is less clear. To examine this we estimate 
fitting limits from the single-band data, indicated by the vertical 
green lines in each panel of Fig.[T3] These lines are defined as being 
0.5 magnitude fainter than the peak in the magnitude histogram for 
Mode_S 1 detections. The blue number in the right column shows 
the number of galaxies brighter than this limit that returned valid 
results in Mode_S 1 fitting, the red number shows the coiTespond- 
ing number using Mode3I fitting. As one can see, the multi-band 
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Figure 13. Magnitude fitting accuracy as a function of simulated magnitude. The blue and red dashed lines show the median scatter of the distribution as a 
function of magnitude for single- and multi-band fits for a common sample of objects for direct comparison. Vertical lines show magnitude limits used in the 
analysis and to derive the objects numbers shown in the different panels. These object counts show the sizes of different samples for which useful values can 
be derived. Detailed description and discussion of these numbers and this figure can be found in the text in Section|5]2] starting on page 12 II 
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fitting returns useful results for appreciably more galaxies, even in 
this bright sample of objects. This illustrates that the increase in 
overall sample size is not only due to adding fainter galaxies, but 
also due to deriving valid fitting results more often for bright galax- 
ies. 

To investigate this topic further, we print two fractions in each 
panel in green. The top number shows the fraction of galaxies simu- 
lated brighter than the single-band fitting limit (vertical green lines) 
for which we obtained valid fits. However, this fraction includes the 
issue of single-band versus multi-band detection completeness and 
is therefore not a direct comparison of the robustness of the fitting 
techniques. In order to take out the detection issue, we determine 
the fraction of objects detected in each band (and above the fit- 
ting limits) for which we obtained valid fitting results. These are 
given by the lower green percentage in each panel. These fractions 
are higher and correspond more closely to the numbers given in 
Table [3] for the full sample. They are also similar between single- 
and multi-band techniques, indicating that for bright samples, at 
least in the case of single-Sersic fitting, the multi-band technique 
does not make large differences to the number of objects which 
can be successfully fit. In some bands this success fraction is actu- 
ally higher for single-band fitting than for multi-band fitting. This 
probably reflects the fact that the cleaning process for multi-band is 
more stringent than that for single band: proximity to a single con- 
straint in any band results in rejection of that multi-band fit. Note 
that even though the individual band success rates are comparable 
between single and multi-band methods, when one desires a sample 
with a complete set of measurements across all bands, the single- 
band method suffers dramatically. This is clearly demonstrated in 
the final three rows of Table|3] It is also important to remember, as 
shown throughout this paper, that the accuracy of the recovered pa- 
rameters, particularly in lower S/N bands, is significantly improved 
by adopting the multi-band technique. 

One effective approach to improving the fractions of detec- 
tions and valid fits is repeat attempts. In the case of detection, if one 
knows an object is present in one band, one can attempt to detect it 
in other bands by varying detection thresholds and deblending pa- 
rameters. For fitting, whenever a fit ends up on a constraint, one can 
restart the fit with slightly different initial conditions or masks, hop- 
ing that one of these repeats will produce acceptable values. Both 
of these approaches are adopted by KI2, with the result that their 
fractions of valid fits are significantly higher than ours. However, 
the fact that restarting GALFIT with slightly different parameters 
can produce different final results is concerning, and it is uncer- 
tain whether parameter values obtained after several attempts are 
reliable. In our case we take constraint violations as a signal that 
reliable results cannot be obtained for the object in question. It is 
not yet clear which approach is most appropriate, although we aim 
to investigate this question further in future work. 

In Figs.[T4]to[T6]we show the equivalent analysis to that pre- 
sented above, but now considering the recoverability of effective 
radii in the simulated data. This complements the similar analy- 
sis that was performed using the results of fitting real data in Sec- 
tion |4.3l As was seen in Section|4l our multi-band fitting technique 
dramatically improves the fitting accuracy for galaxy sizes. This 
is true in all bands (but especially those with low S/N), and at all 
magnitudes. In Figs. [14] to [16] we qualitatively confirm the result 
that we found for real galaxies. However, using the simulations we 
can make more quantitative statements. Fig.[T5] in particular, shows 
that galaxy half-light radius can be measured with much higher pre- 
cision in case of Mode_M compared to Mode.S 1 . Values from the 
top panel in Fig. [15] are given in Table [5] Figure [16] shows that 
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Figure 14. Size recoverability for simulated data (top panel: Mode.S 1 and 
Mode_M, bottom panel: Mode_S2 and Mode_M). Fitting results ai'e im- 
proved when using multi-band fitting. The increased en'or bars in case of 
the entire multi-band sample (in orange), largely come from adding in much 
fainter galaxies, see Fig. I16l and discussion in the text. Only a very small 
systematic offset is visible, even when including the faintest galaxies. 



Table 5. Mean and scatter of the effective radius offsets for Mode.S 1 and 
Mode_M for a sample of objects with complete Mode.S 1 measurements, as 
shown in the top panel of Figure [Tsl 



Band 


mean ztcr Mode.S 1 


mean ±cr Mode_M 


u 


1.06±0.20 


1.00±0.06 


g 


0.99±0.05 


1.00±0.04 


r 


1.00±0.03 


1.00±0.03 


i 


1.00±0.05 


1.00±0.03 


z 


1.04±0.15 


1.01 ±0.04 


Y 


1.04±0.12 


1.02±0.05 


J 


1.04±0.17 


1.03±0.08 


H 


1.03±0.13 


1.05±0.10 


K 


1.03±0.18 


1.02±0.14 



multi-band fitting can measure reasonably accurate sizes even for 
very faint objects, even in bands where the object may be below the 
single-band detection limit. 

Figs.|17|to|19|repeat the same analysis for Sersic indices. Sim- 
ilar to half-light radii, the improvement in fitting quality is much 
clearer than for magnitudes. Be aware that the error bars in all fig- 
ures are relative errors, but they combine values from galaxies with 
different Sersic indices. A galaxy with n = 4 is much harder to fit 
than a galaxy with n = 1, mostly due to uncertainties in the sky 
estimation (H07), so this way of plotting reduces this effect. Val- 
ues from the top panel in Fig. [18] are shown in Table [6] showing 
that the measurements of Sersic index are significantly improved 
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Figure 16. Simulation results showing the recoverability of effective radius for each band as a function of simulated magnitude. Single-band fitting results are 
shown on the left, multi-band results are on the right. Points and lines are the same as for Fig. 1131 The multi-band technique, while only slightly improving on 
magnitude fitting, greatly improves the fitting accuracy for galaxy size. 
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Figure 15. Size recoverability for simulated data, similar to Fig. 1141 but 
only for the subset of galaxies where all Mode.Sl fits (top) or griYHK 
band fits (bottom) return valid results. It is evident that multi-band fitting, 
when comparing the same sample of galaxies, decreases the scatter and thus 
hugely improves the recoverability of galaxy sizes. 
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Figure 17. Sersic index recoverability for the individual bands (top panel: 
Mode.Sl and Mode_M, bottom panel: Mode_S2 and Mode_M). As for ef- 
fective radius, multi-band fitting reduces the scatter on recovered parame- 
ters substantially, without introducing significant systematic effects. 



Table 6. Mean and scatter of Sersic index offsets for Mode.S 1 and Mode_M 
for a sample of objects with complete Mode.S 1 measurements, as shown in 
the top panel of Figure [Tsl 



Band 


mean ztcr Mode.Sl 


mean ztcr Mode_M 


u 


1.27±0.63 


1.06±0.25 


g 


1.02±0.16 


1.03±0.12 


r 


l.OOitO.ll 


1.02±0.08 


i 


1.01±0.13 


1.02±0.10 


z 


1.06±0.36 


1.04±0.12 


Y 


1.08±0.26 


1.04±0.13 


J 


1.02±0.32 


1.05±0.19 


H 


1.02±0.26 


1.08±0.25 


K 


1.06±0.34 


1.09±0.36 



when multi-band fitting is used. In Fig.[T9]we show the effect as 
a function of simulated magnitude. While small systematic effects 
are visible for fainter galaxies, the improvement in fitting quality 
when using multi-band fitting is evident, especially in the lower 
S/N bands. 

We have shown that using multiple band images simultane- 
ously, with structural parameters that are required to vary smoothly 
as a function of wavelength, can dramatically improve both the ro- 
bustness of fits to faint galaxies and the accuracy of recovered pa- 
rameters for objects of all brightnesses. The examples in Figs. [Tj 
and |9] indicate that the recovered wavelength dependencies do re- 
flect the intrinsic variations of the structural parameters between 



wavebands. We now examine the meaning we can ascribe to the 
measured Chebyshev coefficients more closely. 

The simulated galaxies were created such that the simulated 
values for and n follow second order polynomials, e.g. making 
a second order polynomial the perfect fit to the simulated data. We 
have intentionally allowed more freedom, using third-order poly- 
nomials in both Tc and n, when fitting these data with the multi- 
band technique, in order to test how well we can recover the cor- 
rect trends. In Fig.|20]we show histograms of the Chebyshev coef- 
ficients for roi qi (black), q2 (green) and ga (red). There are sev- 
eral things to mention about this figure. Firstly, the histogram for 
gi is clearly offset from 0. As q\ describes the linear terms of the 
Chebyshev polynomials, this means that statistically a non-flat lin- 
ear trend (in this case a negative trend) of with wavelength is vis- 
ible. Galaxies in redder bands generally appear smaller. This effect 
has already been shown in Fig.|5]for real galaxies. As the simula- 
tions were created following this catalogue, the same trend is ap- 
parent in the simulated data, but it is pleasing that we can see this 
directly in the Chebyshev coefficients. Secondly, the histograms for 
both g2 and ga are not significantly offset from 0, so no systematic 
trend is visible in the population. This does not, however, mean 
that individual objects do not significantly show second order poly- 
nomials. Thirdly, and most importantly, the width - and hence the 
importance - of the third order term ga is significantly smaller than 
for g2. Ideally, we would want ga = for all galaxies, but given 
the noise properties of the image, this was not expected. 

A more sensitive way to show this is via the histogram of 
Iog(|g3/g2|) in Fig. [21] Plotted this way the second order term 
'dominates' the shape of the polynomial compared to the third or- 
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Figure 21. Histogram of log([g3/g2|) for galaxies with magsim < 19-8 
This figure compares the importance of the second-order and third-order 
terms of the Chebyshev polynomials for Tc- As the simulations were done 
to second order specifications, one would expect to see that the second order 
term dominates in most cases. The histogram is clearly offset to the left, 
confirmed by the values printed in the figure, which give the numbers either 
side of zero. 



u g r 1 z Y J 



H 



K 



Figure 18. This is the same figure as the top panel of Fig. fTTl but now for the 
subset of galaxies for which all (top) or most (bottom) single-band fits return 
a vaHd result. Again, the reduction in scatter for Mode_M is significant. 



1.0 



0.8 



0.6 



0.4 



0.2 



0.0 



mean 
-0.61 

0.09 
0.01 




width 
1.59 




-4 



q|1,2,3],r( 



Figure 20. Histogram of the Chebyshev coefficients describing the wave- 
length dependence of r'c : gi (black), q2 (green) and 93 (red), obtained from 
G ALFfTM multi-band fits to our simulated dataset. The mean and rms scatter 
of each distribution are indicated by vertical solid and dashed lines, respec- 
tively, and the values printed in the comers of the figure. Histograms have 
been normalized to unit maximum. As expected, the distribution of 93 in 
our simulated sample is narrower than the distribution of 92 ■ 



der term in all galaxies that show log(|g3 /q2 1) < 0, or equivalently 
I93I < I §2 1 ■ This is important in order to justify whether the polyno- 
mial order fitted to the real data is the correct order to use - e.g., if 
one fits a third order polynomial and finds that |g3 1 < |g2 1 for most 
galaxies, it is a justifiable assumption that second order polynomi- 
als are appropriate for most objects. In previous tests this is what 
we have found, hence our decision to fit real galaxies with second 
order polynomials and creating simulations to this specification. 



6 COMPUTATIONAL CONSIDERATIONS 

Other than fitting accuracy, there are also other issues to be con- 
sidered. In this section, we will discuss some of these issues and 
compare the single and multi-band fitting procedures. 



6.1 Fitting time 

Potentially, multi-band fitting may require less computing time to 
produce useful fits for large samples of galaxies, when compared 
with the single-band approach, as in principle we have fewer free 
parameters being constrained by the same dataset. As we have fit 
samples using both methods, we can directly make this comparison. 
Here we consider only the fitting time, i.e., the processing time 
required to process the galaxy with GALFITM. We ignore the time 
it takes for the fit to be initially set-up, as this is similar for both 
methods. GALAPAGOS does most jobs on a band-by-band basis, so 
the required setup time should generally be the same. However, for 
the multi-band method not all steps in the code are repeated for all 
bands; e.g. deblending/masking decisions are common for all bands 
and only decided once, whereas Mode_S 1 makes these decisions on 
each band independently. Some further optimization of this part of 
GALAPAGOS-2 is also possible in the future. 

When comparing the pure fitting time for the simulated ob- 
jects, we find that the multi-band fits are in fact slower than the 
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Figure 22. CPU time histograms. Red liistograms represent times for 
Mode_M fitting, black fiistograms give tlie sum of tlie times for the nine 
Mode.Sl fits. This figure only contains values from objects where all nine 
Mode_S2 and Mode_M fits exist, to allow a fair compaiison between the 
methods. Vertical lines show the mean and rms scatter for both histograms. 
The top panel shows the total fitting time as returned by GALFITM. Multi- 
band fits take longer than single-band fits. The main reason for this is seen 
in the second panel, every individual iteration takes substantially longer in 
Mode_M compared to Mode_S2. This effect is nearly, but not entirely, can- 
celed out by the smaller total number of iterations needed to derive all 9 
band fits (lower panel). 



nine individual band fits combined. The top panel of Fig.l22lshows 
a histogram of the overall fitting time for a small subset of objects 
for which we have all nine individual band fits, and for which such a 
comparison is possible. The fitting time will depend on how many 
secondary objects with free parameters, i.e. those fainter than the 
primary, are included in the fits. The number of secondaries de- 
pends upon the deblending, and is greater for multi-band detec- 
tion. To compare fitting times fairly, we therefore plot results for 
Mode_S2 and Mode3I. While the combined single-band fits take 
around 74 seconds/object, the multi-band fits take about 106 sec- 
onds/object. 

The most time consuming part of the GALHTM algorithm is 
producing the model derivatives, which describe how the model im- 
ages would change for small variations in each of the parameters. 
The derivatives must be calculated for every pixel and convolved 
with the PSF, which is an expensive operation. The Levenberg- 
Marquardt algorithm requires these derivative images, one set for 
each band being used in the fit, to be calculated with respect to 
each parameter, at every iteration step. We would therefore naively 
assume that the fitting time per iteration should be roughly propor- 
tional to the number of free parameters. 

For instance, a single-band, single-Sersic fit has 7 free param- 
eters (magnitude, a;-centre, j/-centre, half-light radius, Sersic index, 
axis ratio and position angle), and hence 7 derivative images must 
be calculated by GALFITM each iteration. However, this is done in- 
dependently for the 9 input images. The multi-band fits caiTied out 
for this work use 19 free parameters (9 for magnitude, 1 x-centre, 
1 y-centre, 3 for half-light radius, 3 for Sersic index, 1 position 
angle and 1 axis ratio), and derivatives must be calculated for 9 
bands. Hence, 19 x 9 = 171 derivatives must to be calculated for 
each iteration. By this estimate, each multi-band iteration should 
take 171/7 ~ 24 times longer than each single-band iteration (this 
factor should be the same for any number of secondary objects in- 
cluded in the fit). In reality, from the middle panel of Fig.|22l we 
see that the difference is somewhat less: a factor of ~ 13. This is 
probably because the derivatives of higher order Chebyshev param- 
eters are less expensive to calculate than the zeroth-order parame- 
ters and, similarly, the standard parameters in single-band fitting. 

When one accounts for the fact that the single-band fit must be 
run for each of the bands, the multi-band method gains by a factor 
of 9. If the number of iterations required for each method were 
the same, the multi-band method would take 13/9 times as long 
and hence would be ~ 40 per cent slower. Indeed the number of 
iterations required for each fit are very similar, with the single-band 
approach requiring about 9 times more iterations, as the fit must be 
performed for each band, as shown in the lower panel of Fig. 122 1 . 
For these bright galaxies, with complete sets of single-band fits, 
the multi-band method therefore takes roughly 40 per cent more 
fitting time compare to the Mode_S2 single-band fits. Note that, 
for single-band detection, the number of secondaries that are fit is 
reduced, and hence the combined time of the Mode.Sl fits take 
only ~ 39 seconds compared to ~ 81 seconds for the same sample 
of galaxies in Mode.M. The multi-band Mode_M method therefore 
takes roughly twice the time of Mode.S 1 . 

The number of successfully fit objects is much higher in the 
case of multi-band fitting (see Table[3}. From the above arguments, 
this should result in a much larger number of galaxies with a com- 
plete set of multi-band measurements per unit processing time. 
However, the above speed comparisons are limited to the cases 
where we obtain successful fits in all nine single-band fits and the 
multi-band fit. For these, the number of iterations are similar, but 
for the majority of galaxies, with lower signal-to-noise, the single- 
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band fits (many of which fail) appear to require significantly fewer 
iterations than the multi-band fits. To get an overall impression of 
the speed of single- versus multi-band, we divide the number of 
galaxies with a complete set of multi-band measurements by the 
total time taken to fit the entire sample, to produce a 'complete suc- 
cess' rate. For Mode_S2 this is ~ 2.4 'complete successes' per hour 
(if one requires all ugrizY J H K band values; 23.5 requiring just 
the griYHK bands), while for ModeJVI it is 29.3 'complete suc- 
cesses' per hour (full ugrizY JHK). In single-band fitting, much 
more CPU time is spent on objects that are later discarded as un- 
reliable by the cleaning process. So in terms of useful fits per unit 
time, multi-band fitting compares very favourably. 

6.2 Storage requirements 

Another practical aspect to take into account, especially for large 
datasets, is the disk space used. If a dataset uses less disk space it is 
easier to handle and less space will be needed for big surveys. For 
the sample of 29205 real objects, we use: 

• 51 GB for the Mode.M fitting, 

• 78 GB for the Mode.Sl fitting. The difference is largely down 
to having SExtractor files for each band instead of one set for the 
whole process. In comparison, the single-band fitting does save 
some disk space because not as many galaxy postage stamps have 
to be made for the individual bands as less objects are detected. 

• In Mode_S2, cutting all postage stamps when running the ob- 
ject detection on a multi-wavelength image increases the disk space 
required to a total of 100GB, so a reduction of disk space ^ 50 per 
cent is achieved when changing from single-band to multi-band fit- 
ting. 

Multi-wavelength fitting has the smallest space-to-result ratio, es- 
pecially when taking into account that one gets usable results for 
a much bigger sample. We get valid results (fits that did not run 
into fitting constraints) for 15566 objects, compared to 209 galaxies 
with valid fits in u-band, and 3276 (ff-band), 5802 (r-band), 6804 
(i-band), 2117 (z-band), 3912 (F-band), 2530 (J-band), 5432 {H- 
band) and 5190 (_ftr-band), respectively. The overlap of all these 
samples is a mere 87 galaxies (or 992 when only griYHK are 
taken into account). 

For the set of simulated galaxies, we find a similar result: For 
a total of 42229 simulated objects, the disk space needed is: 

• 64 GB for Mode_M fitting, 

• 73 GB for Mode.Sl fitting, 

• 1 14GB for Mode_S2 fitting. 

Object numbers with useful values here are: 714 (u-band), 6022 
(S-band), 10362 (r-band), 12326 (i-band), 4139 (z-band), 9627 
(y-band), 7965 (J-band), 12593 (//-band) and 12602 (/f-band) 
with an overlap of 305 objects (2560 objects). This compares to 
23632 objects for which we get a valid result in all bands when 
using multi-band fitting. Overall, our multi-band fitting procedure 
can help to save nearly 50 per cent of disk space required. 



7 COLOUR-MAGNITUDE DIAGRAM FROM 
MULTI-BAND FITTING 

The multi-wavelength Mode_M technique that we have presented 
and tested, both above and in accompanying MegaMorph papers, 
has subsequently been applied to a larger sample of galaxies from 



the GAMA survey. Here we use this sample to present prelimi- 
nary colour-magnitude and size-magnitude diagrams. This is done 
for two purposes: as a consistency check of the fitting results, and 
to demonstrate where the multi-band technique offers scope for 
quickly improving science results. 

Single-band single-Sersic fits, in ugrizY J H K , have been 
performed by K12 for all spectroscopic targets in the GAMA sam- 
ple, using their own code, SIGMA, to run GALFIT3. However, a 
detailed comparison between our single-band fits and those con- 
ducted in K12 is beyond the scope of this paper We therefore de- 
cided that comparing their single-band fitting results to our multi- 
band results could potentially lead to wrong conclusions regard- 
ing single vs multi-band fitting, due to the use of different wrapper 
scripts. However, in earlier tests, we have established that both pro- 
cedures, SIGMA (K 12) and our own method based on GALAPAGOS 
tearden et alj201 j) generally return very similar results when used 
on single-band data, and comparable values when used on multi- 
band data in our setup. In order to show the strength of our code, 
we avoid using any additional information about the objects other 
than GALAPAGOS fitting results and spectroscopic redshifts. The 
latter are provided by GAMA, and used only to calculate physi- 
cally meaningful parameters such as absolute magnitudes and sizes, 
for both our single and multi-band fitting results. For these trans- 
formations, a cosmological model of Q./^ = 0.7, f7m = 0.3 and 
//o — 70 hji) km s^^ Mpc~^ is assumed. 

In Section |4] we have carried out both single and multi-band 
fits on a small subarea of the GAMA survey field. Compared to the 
GALAPAGOS version used and tested above, we use a minimally 
changed version in this section. The main change was the intro- 
duction of a feature that allows the code to only target certain ob- 
jects instead of targeting all detections in a survey. For this chapter 
we applied (using Mode_M) multi-band version of GALAPAGOS-2 
and GALFITM to all objects that have both full wavelength cov- 
erage in the GA MA 9-hour (G09) fi eld and redshifts provided by 
the GAMA team jPriveret al]|201lL Liske, in prep). Additionally, 
we include all objects within 66 arcsec (~ 200 pixels) to one of 
these GAMA objects and brighter than said object, as GALAPAGOS 
works on objects in brightness order from bright to fainll"1. A good 
fraction of these neighbouring objects are bright stars. Given the 
brightness of the objects targeted, this slightly changed code should 
not influence fitting results significantly but allows the code to only 
target 144261 objects instead of the total 1401969 detections in the 
field, reducing the sample size by ~ 90 per cent and saving a huge 
amount of fitting time. 

Of these 144261 objects, we have spectroscopic redshifts for 
43617 objects. GALAPAGOS returned a result for 43572 (99.9 per 
cent), and 33696 (77.3 per cent) have useful fits and a full set of 
multi-wavelength parameters is provided. The fraction of galaxies 
with valid fits is significantly higher in this sample of galaxies com- 
pared to the samples used earlier in this paper (e.g. the left part of 
Table |2j, mostly because for this chapter bright galaxies - that are 
potentially easier to fit - have been primarily targeted, whereas in 
the previous chapters, we had targeted all detected objects. Also, as 
stars have no redshifts assigned, they are removed when deriving 
these numbers, boosting the success rate. 

The resulting catalogue of this procedure is used in both Vul- 



During the development of the code, it proved useful for fitting stabil- 
ity if results for brighter neighbouring objects are already known at the 
ti me they have to be u sed as neighbours. This issue is discussed in detail 
in lBardenet"ZI)20I2l) . 
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Figure 23. Colour-magnitude diagram (y-axis shows purely observed colours). Top-left: when using aperture photometry (GAMA SExtractorMAG_AUTO 
from catalogue ApMatchedCatv03); Top-right: when using single band fitting (Mode.Sl). colour-coding using single-band values; Bottom-left: Mode.Sl 
fitting results, colour coding using Mode_M riz values. Bottom-light: using only Mode_M multi-band fitting, colour coded by z-band Sersic index from multi- 
band fitting (blue=low (< 0.5), red=high(> 6)). Uz is assumed to be the best separator and we argue in this paper for multi-band to be a more reliable value 
than single-band Sersic indices. A separation between 'early-type' (high-n) and 'late-type' (low-n) galaxies becomes much clearer when ModeJVI values 
are used for colour coding. Please be reminded that the number of points in this analysis is mostly defined by the number of redshifts available, and not by the 
number of successful fits. As a consequence, single and multi-band figures contain a very similar number of points. 



cani et al (2012, in prep) and this section, in which we briefly ex- 
amine the colour-magnitude diagram of galaxies. 

In Fig.|23l we show g-i vs. Mi Colour-Magnitude Diagrams 
(CMDs) as they would be recovered by the different techniques, 
colour-coded by Sersic index n^. The top-left panel shows values 
from aperture photometry, for which we use matched-ape rture pho- 
tometry obtained in a manner similar to lHill et^ I( |201 l|) and to be 
described in Liske et al. (in prep). The points in this panel have 
not been colour coded as this information would have to be pro- 
vided by a second code, aperture photometry itself does not give an 
estimate for galaxy 'morphology'. The top-right figure shows the 
CMD as given entirely by Mode_S 1 fitting, e.g. the colour coding 
according to as provided by single-band fitting. High-n and low- 
n galaxies show mixed colours in this figure, a separation, although 
indicated, is weak. The bottom-left figure shows the same single- 
band values, but colour-coded using multi-band n^. In this figure, 
the separation between low-n and high-n galaxies is better. The 
reader should keep in mind that Sersic index is not a robust value 
to be used for classifying galaxies in a morphological sense, espe- 
cially early-type galaxies are known to show a wide range of Sersic 
indices, e.g., they are following a magnitude-n relation (see e.g. 



lGrahanill20T ll). Additionally, Sersic indices are sensitive to mea- 
surement errors, for example, a faint AGN or bright star cluster in 
the centre of a galaxy would boost n to show a high value instead of 
the intrinsic value of the underlying host galaxy. We use n^ only as 
a very crude estimator in this analysis and would strongly discour- 
age users to use it for more than that. Especially, it cannot replace 
a visual classification of morphology. 

In the bottom-right panel of the figure we show the CMD as 
recovered entirely by multi-band fitting. The division between low- 
n and high-n galaxies becomes much stronger, and the scatter in 
the figure gets reduced, with high-n galaxies showing significantly 
redder colours than low-n galaxies. This is expected by simple 
morphology-stellar-population arguments. It should be stated here 
that galaxies are plotted in random order, not by increasing Sersic 
index, the lack of blue points on the red-sequence and red points in 
the blue cloud, even when using n^ as colour index, is real and not 
due to plotting procedure. 

Please note that the number of points in the single-band panels 
(727 objects) is only slightly smaller than the one in the multi-band 
panel (795 objects). This is true as long as no u- or z-band data is 
used and only bright galaxies which were targeted for spectroscopic 
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redshifts are considered (m^ < 19.8). Given that this completeness 
effect has already been discussed in detail above, we decided to not 
use u- or z-band data here, although it might show a better separa- 
tion of the red sequence and the blue cloud. Instead we concentrate 
on the fitting accuracy itself. 

Given the overall shape of this figure, it should be empha- 
sized that this figure shows purely observed colours, e.g., no k- 
corrections or dust-corrections have been applied. The main reason 
for this is that we want to show the difference between the tech- 
niques, not sample sizes. Robust k-corrections would require mag- 
nitude measurements in many bands which are automatically pro- 
vided by multi-band fitting. Single-band fitting, however, can pro- 
duce results in some bands, not in the others, providing reliable k- 
correction only for a subset of galaxies. We will apply k-correction 
to all multi-band fits in Fig. 1241 

In Fig.|24]we show the CMD for the 28737 GAMA galax- 
ies in the G09 field for which we were able to recover useful fit- 
ting values using multi-band fitting procedures (Mode_M). In con- 
trast to Fig. |23] values in this figure have been k-corrected to show 
restframe values, using KCORRECT dBlanton & Roweisir2007h . ver- 
sion 4_2. This is possible here because the multi-band fitting pro- 
cedure returns magnitudes for all ugrizYJHK bands, whereas 
any single-band routine could have some magnitude values miss- 
ing in certain bands, making k-correction both harder and less ac- 
curate. The middle panel shows contours of the same distribution 
for all objects. An indication for a bimodality and galaxy separa- 
tion can be seen. The right panel shows the same, but galaxies have 
been split up into broad -value bins, each with a large number 
of galaxies (7000-12000 galaxies), in order to show the separation 
between the different Sersic index populations more clearly. The 
bimodality becomes very visible here, although the red sequence is 
not as tight as one would hope when a clean classification scheme 
is used instead of values. 

A few remarks should be made on these findings: 

• The solid line in the rightmost panel of Fig. |24] indicates not 
a fit to a red-seque nce, but a separator, a fit to the 'green valley', 
as determined by iGavazzi et alj ( 1201 Ol) from a sample of 4100 
galaxies using visual classification. Gavazzi et alj ( 1201 Oh use SDSS 
DR7 Petrosian magnitudes for their analysis. Due to the nature 
of these - the r-band alone defines the size of the aperture in 
which the magnitude is measured - we would expect a widening 
of our red sequence due to the ability of our method to vary galaxy 
size with wavelength, hence measuring a slightly di fferent magni- 
tude a t other wavelengths than the ones provided bv lGavazzi et alj 
( I2OI O"). We argue that our magnitudes more closely resemble the 
true magnitudes at these wavelengths in comparison with SDSS 
Petrosian magnitudes which should be bulge-dominated given that 
the aperture size is defined in a reddish (r) band (please also see 
iGraham et al.l2005h . Unfortunately for this comparison, this bulge- 
domination of the magnitudes measured helps to get a tight colour- 
magnitude relation and a better separation from the blue cloud. 

• It is very interesting to see that the contours for galaxies with 
intermediate (1.5 < n < 3; shown in green) seem to have a 
peak on the red sequence and a long 'tail' towards the blue cloud. 
This is exactly what would be expected for a transition population. 

• Our galaxy 'separator' Uz should only give a very rough esti- 
mation of the galaxy 'morphology' as it is easily biased and boosted 
by central light sources such as faint AGN and/or bright star clus- 
ters. Hence the high-n sample will include some of these objects, 
although they would visually be classified as disks, suggesting a 
place in or closer to the blue cloud. A way to test this hypothesis is 
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Figure 25. The magnitude-size relation. Points are colour coded by riz 
matching the previous s (blue=low (< 0.5), red=high(> 6)). As discussed 
in the text, there is a visible population of very bright galaxies with in- 
creased riz and very small sizes, visible as a cloud of red points at the 
lower-left edge of the cloud of blue points. 



to look at the magnitude- size relation of all objects. When a single 
Sersic profile is used to fit an object which contains both a disk and 
a point source, one would expect to fit both a high Sersic index and 
a small size, as well as a bright magnitude as this would include the 
flux of the nuclear source. When plotting the magnitude-size rela- 
tion of all the objects in our sample (see Fig. |25t , we can see that 
indeed many of these high-n galaxies are the brightest and small- 
est galaxies, seen as a purple/red cloud of points, in what would 
traditionally be called the 'blue cloud' of disk-dominated galaxies. 

• Even using perfect classifications, one would always find 
both blue ellipticals in the blue cloud and disks that show red 
colours given their dust properties and possibly seen edge-on (e.g. 
[Rowlands et"ai]|2012l ; IWolf et alj|2009l) . An entirely clean separa- 
tion is impossible to achieve. 

• Please again bear in mind that no other correction (e.g., for 
dust content) than k-correction has been applied to any of the values 
in Fig. [24] This preliminary analysis is intended only as a demon- 
stration of the advantages of multi-band fitting. 

Our simple analysis, using no other information than redshifts 
and colours as directly measured by our multi-band fitting code, 
clearly shows the known separation between the red sequence and 
the blue cloud of galaxies, even when using a crude and inaccurate 
proxy for 'galaxy morphology' such as Sersic index Uz- Although 
a separation of the two galaxy populations by using Sersic index as 
a 'classifier' will never be ideal, it proves a powerful statistical tool 
in this example. 

We assume that most of the outliers in Fig. [24] are due to an 
unflagged eiTor/uncertainty in our fitting pipeline and our quick and 
rough analysis. Especially, we are not (yet) running bulge-disk de- 
composition on our data, and fitting single Sersic profiles to B/D 
composite galaxies is known to be non-ideal. We assume that the 
fitting process went wrong in these examples in a way that we 
were unable to detect in a simple automated fashion. Statistically, 
these ~ 100-200 outliers become insignificant in a sample of 28737 
galaxies. 
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Figure 24. Colour magnitude diagram (botli axes restframe and k-corrected). Left: all GAMA G09 galaxies with redshifts; middle: contours for all galaxies in 
the left figure; right: contours for samples split by Sersic index. Points are colour coded by z-band Sersic index (blue=low (< 0.5), red=high(> 6)). A clear 
separation between low-n and high-n galaxies is visible (please note that the points are plotted in random order, not by increasing Sersic index), especially 
when plotting three n bins (n < 1.5, 1.5 < n < 3, 4 < n) as contours (right figure). Even in the contours for all galaxies (Center), the shape of the contours 
indicate 2 overlapping sequences, which is what one ex pects. The solid Une in the light figure indicates a separator between the red sequence and the blue 
cloud as empirically determined bv lOavazzi et al 



8 SUMMARY 

In this paper, we have presented new versions of both GALFIT 
(which we call GALHTM) and GALAPAGOS (GALAPAGOS-2) to au- 
tomate the process of fitting two-dimensional, single-Sersic mod- 
els to large samples of galaxies. Our novel approach uses multi- 
ple images taken at different wavelengths simultaneously, instead 
of single-band data, as used by the current published versions of 
these codes. We have tested our multi-band modifications exten- 
sively, both on real and simulated data and compared their perfor- 
mance to the single-band versions of the codes. For real data, we 
have selected 9-band ( uqrizYJHK) data from the GAMA sur- 
vey jPriver et al .11201 1 ). Simulations were created following H07 
to largely reproduce the parameter distributions of the real data. 
This creates a dataset as realistic as necessary, while taking out the 
issues of image confusion by nearby stars and internal galaxy struc- 
ture, as well as providing a truth value for comparison against those 
recovered recovered by the fitting process. 

From both real and simulated data, we can draw very similar 
conclusions: 

• Multi-band fitting significantly improves the recovery of 
galaxy magnitudes, at least in the lower S/N bands (particularly 
u and 2-band in this analysis, see e.g. Fig.llOland Table|4l(. The im- 
provement for both sizes (see e.g. Fig.[T4]and Table|5) and Sersic 
indices (see e.g. Fig. [17] and Table [SJl is significant over the entire 
wavelength range tested in this paper. The multi-band approach re- 
turns more accurate results than treating the bands independently, 
and is able to make useful measurements for very faint galaxies. In 
our analysis, multi-band fitting is able to derive reliable values for 
~ 75 times more galaxies in SDSS u-band data than the single- 
band approach. 

• The multi-band version of our code is more stable, in terms 
of providing fitting results for objects where single-band fits fail 
to return a result. This is because the single-band fits occasionally 
return results which are affected by the constraints that are provided 
by GALAPAGOS and used by GALFITM during the fitting process. 
GALFITM more frequently returns a valid results when multi-band 



data is used. This is always the case for real galaxy images, while in 
our simulations this is true when fits in several bands are required, 
e.g. for SED modelling. 

• Physically meaningful variations of galaxy parameters with 
observed wavelength are used in the fitting method presented in 
this paper. They enable accurate wavelength interpolation, useful 
to derive restframe galaxy parameters, and provide the user with 
a full set of parameters at all wavelengths. Single-band fits on the 
other hand are more likely to fail in at least individual bands, while 
returning a result in other wavelengths, thus leaving a less complete 
sample from which to interpolate values. 

• While running single-band fits returns more or less equivalent 
results for bright galaxies, multi-band fitting provides a more ho- 
mogenous dataset as decisions about object detection, deblending 
and masking are necessarily the same. Single-band fits can vary 
these decisions between different bands. 

• In addition to providing results for a higher fraction of galax- 
ies, the multi-band code can be run on fainter objects, as the object 
detection can be sensibly carried out on a co-added image. The 
number of galaxies that can be used for scientific purposes is in- 
creased by a factor of ~ 2 if only single-band values are needed 
up to ~ 16 if multi-band data {griYHK) values are required ~ 
a factor of ^ 180 is achieved when a full set of parameters in 
ugrizYJHK data are needed. 

Over all, we conclude that using multi-band fitting instead of 
single-band fitting improves both quality and quantity of galaxy fits 
on survey data. It therefore provides an excellent tool for modem 
surveys and enables a step forward in many areas of astronomy, es- 
pecially where sample size and magnitude constraints limit current 
work. 



9 PROSPECTS 

In this paper, we have exclusively tested GALFITM and 
GALAPAGOS-2 on single-Sersic profiles. Further plans most im- 
portantly include implementation of bulge-disk decomposition into 
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GALAPAGOS-2. As some of the conclusions in this paper might 
change due to both the additional flexibility and the more challeng- 
ing fits even on bright objects, we will discuss these changes in 
detail in the follow-up paper. The intention is to separate galaxy 
bulges and disks more accurately than previously possible by using 
the information provided by the different colours of the respective 
components. 

The colour gradients (variations of n and re with wave- 
length) seen in this paper come from a combination of both mix- 
ing the bulge and disk stellar populations when fi tting them with a 
single-Sersic profile and the effects of dust (e.g.. lMollenhoff et alj 
l200d|Pastrav et al.ll2012h . Fully understanding the observed trends 
will require careful consideration of both multi-wavelength bulge- 
disk deco mpositions and spatially-resolved models of dust attenu- 
ation (e.g.. |Popescu et al][2011 ). constrained by spacebome UV & 
FIR/submm photometry." 

Another issue, which we have tried to encompass by testing 
with both real and simulated data in this work, but which has other- 
wise been neglected, is the performance of multi-band GALFITM 
(and profile fitting in general) in the presence of morphological 
features, such as bars and spiral arms. We aim to investigate this 
in fu ture by drawing on the data and resources of Galaxy Zoo 
jTintott et al. 2011 ; Masters et al. 201 1) 

Although GALFITM already supports multiple components, an 
adapted version of GALAPAGOS-2 will be needed for this work. 
This version will be presented in a different paper (HauBler et al., 
in prep). It is currently under development and testing and we aim 
for it to be ready for public release in ^^6 months. We will publish 
our codes and tools in a follow-up paper when bulge-disk decom- 
position is successfully incorporated into the code and encourage 
astronomers to use these tools on their data. 
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APPENDIX A: EXAMPLE IMAGES AND FIT RESULTS 

In this appendix, we show some sample images of simulated galax- 
ies, their fitting models (if applicable) and their fitting residuals 
in all 9 bands for both multi-band and single-band fitting. These 
galaxies are not randomly chosen, but rather chosen to show typi- 
cal difficulties and effects that separate single-band and multi-band 
fitting. In particular, the objects shown are the same objects shown 
in Fig.[2]andFig.|9l 
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Figure A2. Example real object, similar to Fig. lAll This example shows a galaxy that has not been detected in single-band M-band. The neighbouring galaxy 
has also not been detected in the single-band g-band, potentially biasing the fit. 
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Figure A6. Example simulated object. From left to right: input image, model & residual multi-band fit, image, model & residual single-band fit. Black images 
in case of single-band fitting indicate where no GALFITM results exist (while ignoring constraints and cleaning of the catalogue). 
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Figure A7. This is an extreme simulated object. It can really only be seen in Y- and J-band data, but thanks to the multi-band capabilities, we can recover 
a full SED and parameter in all bands on this galaxy. In comparison, single-band fitting only returns one value in the J-band fit. Please keep in mind that 
all examples in this Appendix ignore the cleaning of the catalogue, e.g. showing this example does not mean that this objects would end up in the science 
catalogue. 
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